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CHAPTER 


I 


INTRODUCTION 

Interest  in  underwater  sound  transmission  has  led 
to  the  study  of  the  acoustic  properties  of  marine  sedi- 
ments. Extensive  experimental  work  has  established  useful 
empirical  relationships  which  characterize  compressional 
wave  propagation  in  such  materials.  However,  compressional 
wave  parameters  alone  are  inadequate  for  a complete  de- 
scription of  acoustic  interactions  at  the  ocean  floor.  Al- 
most all  natural  sediments  exhibit  sufficient  rigidity  to 
support  transverse  oscillations  excited  by  compressional 
wave  conversion.  The  shear  wave  properties  of  marine  sedi- 

V 

ments,  however,  are  poorly  known.  The  shear  modulus  is 
normally  so  small  that  direct  measurement  of  shear  wave 
properties  has  been  difficult.  Although  scattered  data 
exist,  a general  experimental  characterization  of  the  vis- 
coelastic shear  wave  properties  of  fully  saturated,  uncon- 
solidated media  has  been  lacking.  Recently,  the  develop- 
ment of  new  transducers  has  overcome  the  barriers  to  shear 
wave  detection  in  low  rigidity  materials  (Shirley  and  Hamp- 
ton, 1978;  Shirley,  1978).  Such  transducers  were  utilized 
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in  the  present  study  to  obtain  shear  wave  speed  and  attenu- 
ation in  several  laboratory  sediments. 

A review  in  Chapter  II  of  previous  experimental 
determinations  of  sediment  shear  wave  properties  reveals 
insufficient  data  to  test  the  various  theoretical  models  of 


shear  wave  propagation  reviewed  in  Chapter  III.  Measure- 
ments of  shear  wave  speed  and  attenuation  as  functions  of 
frequency  in  the  same  sediment  sample  are  required  for  a 
decisive  comparison  of  theoretical  predictions.  Various 
physical  properties,  notably  the  permeability,  are  also  re- 
quired. Chapters  IV  through  VII  detail  the  acquisition  of 
such  data  in  the  present  study.  Chapters  VIII  and  IX  pre- 
sent the  results. 

Although  little  experimental  work  has  been  done 
in  surficial  marine  sediments,  interest  in  shear  wave  prop- 
agation in  porous  media  is  not  confined  to  underwater 
acoustics.  The  dynamic  shear  modulus  of  loose  soil  is  a 
determining  factor  in  the  design  of  foundations  subject  to 
vibratory  loads  (Wu,  1971).  Combined  shear  and  compres- 
sional  wave  measurements  can  potentially  provide  the  petrole- 
um industry  with  a direct  indication  of  hydrocarbon  accumu- 
lations (Tatham  and  Stoffa,  1976).  Changes  in  shear  wave 
propagation  caused  by  the  reaction  of  microcracks  to 
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abnormal  stress  are  a possible  means  of  earthquake  predic- 
tion (Nur , 1972).  Theories  of  acoustic  propagation  in  two- 
phase  fluid  and  solid  media  have  been  applied  to  interior 
portions  of  the  earth  wherein  the  pore  fluid  is  molten  rock 
(Walsh,  1968;  Gliko,  1976). 

Numerous  experimental  and  theoretical  studies  of 
shear  wave  propagation  applicable  to  the  above  interests 
have  been  made.  These  cannot  always  be  directly  applied  to 
the  study  of  saturated  sediments.  Partial  saturation, 
large  vibration  amplitudes,  and  large  confining  pressures 
are  all  inappropriate  for  surficial  marine  sediments.  Nev- 
ertheless, important  insights  into  shear  wave  propagation 
can  be  gained  from  such  studies.  These  will  therefore  be 
reviewed  along  with  the  available  sediment  measurements  in 


Chapter  II. 


CHAPTER 


I I 


EXPERIMENTAL  BACKGROUND 

' 

Review  of  Techniques 

Several  experimental  procedures  have  been  devel- 
oped to  measure  the  dynamic  properties  of  porous  media. 
Taken  as  a group,  these  techniques  provide  for  stress  cy- 

. 

cles  from  less  than  a hertz  to  several  megahertz.  However, 
the  individual  methods  have  limited  resolutions  and  fre- 
quency sensitivities.  It  is  therefore  necessary  to  under- 
stand the  limitations  of  the  various  experimental  arrange- 
ments before  comparison  of  their  results  are  undertaken. 

Since  the  properties  of  a grossly  homogeneous  me- 
dium are  desired,  the  size  of  the  available  sample  places 
the  first  restrictions  on  the  measurement  of  its  acoustic 
parameters.  The  concept  of  a porous  substance  as  a single 
entity  is  invalid  when  the  wavelength  approaches  the  pore 
size.  Thus,  the  maximum  frequency  of  interest  is  deter- 
mined by  the  medium  itself.  The  lowest  frequency  available 
to  an  investigator  with  a finite  sample,  either  in  the 
field  or  the  laboratory,  depends  on  the  validity  of 
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techniques  that  account  for  the  impedance  contrasts  at  the 
boundaries . 

Laboratory  Pulse  Transmission 

Direct  laboratory  measurement  of  the  acoustic 
properties  cf  a porous  media  is  easily  obtained  from  the 
propagation  of  an  acoustic  signal  through  a sample  large  in 
comparison  to  the  wavelength.  Generally,  the  transmission 
path  lies  between  a projector  and  one  or  more  receivers. 

The  transit  time  of  a pulse  or  the  phase  change  of  a con- 
tinuous signal  allows  determination  of  the  speed.  A change 
in  signal  amplitude  with  separation  indicates  the  energy 
loss.  Other  arrangements  place  both  projector  and  receiver 
on  one  edge  of  a sample1  and  monitor  the  arrival  times  and 
amplitudes  of  multiple  reflections  from  the  opposite  bound- 
ary. This  provides  a change  in  path  length  for  a single 
fixed  receiver. 

If  broadband  transducers  are  used,  the  variation 
of  amplitude  with  frequency  may  be  monitored  by  either  of 
two  common  methods.  The  ratio  of  * he  Fourier  spectra  of  an 
impulsive  input  recorded  at  two  receivers  gives  the  rela- 
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tive  amplitude  loss  over  a set  propagation  distance. 
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Alternatively,  one  or  more  cycles  of  a sine  wave  may  be 
used  as  the  input.  A narrow  bandpass  filter  centered  on 
the  oscillator  frequency  is  applied  to  the  received  signal 
to  eliminate  frequencies  introduced  by  gating  the  input. 

The  frequency  resolution  is  determined  by  the  filter  set- 
ting and  the  number  of  cycles  present  in  the  pulse.  Ideal- 
ly, perfect  frequency  discrimination  could  be  obtained  from 
a continuous  input.  However,  this  is  normally  impractical 
due  to  interference  from  reflections  present  in  the  wave 
train.  Reflections  can  also  limit  the  application  of  spec- 
tral analysis.  In  any  case,  the  responses  of  the  measuring 
devices  and  the  geometric  divergence  of  the  energy  must  be 
determined  for  a valid  measurement  of  the  attenuation.  De- 
pending on  the  system,  velocities  reported  in  the  litera- 
ture  may  be  group  velocities  or  phase  velocities. 

The  transducers  used  for  laboratory  compression- 
al  wave  measurements  are  normally  piezoelectric  crystals  or 
polarized  ceramics.  The  useful  frequency  output  depends  on 
the  type  of  material,  the  size  of  the  element,  and  the  meth 
od  of  coupling.  Some  experimental  arrangements  use  broad- 
band elements,  others  use  highly  resonant  devices.  Gregory 
(1975)  gives  a short  discussion  of  the  advantages  and  dis- 
advantages of  each  technique.  The  narrowband  signal  offers 


better  time  resolution  and  is  less  subject  to  frequency 
distortion  in  materials  with  high  attenuation.  The  choice 
of  a particular  ceramic  material  is  sometimes  dictated  by 
the  ambient  pressures  encountered  in  the  experiment.  If 
the  element  is  not  isolated  from  pressure  changes,  materi- 
als are  chosen  which  are  relatively  insensitive  to  the 
pressure.  These  normally  are  low  Q devices  (Shirley  and 
Bell,  1978). 

Laboratory  shear  waves  can  be  generated  by  a ce- 
ramic element  polarized  perpendicular  to  the  applied  field 
or  by  quartz  crystals  cut  in  either  the  Y or  AC  planes 
(e.g.,  Simmons,  196U).  Alternatively,  compressional  waves 
may  be  converted  to  shear  waves  at  an  interface  (King, 

1966;  Gregory  and  Podio,  1970).  The  second  method  provides 
a shear  wave  pulse  free  from  the  compressional  wave  precur- 
sor that  is  often  generated  by  shear  elements.  Both  of 
these  techniques  have  allowed  successful  measurements  of 
shear  wave  properties  of  rocks.  However,  such  instruments 
often  do  not  have  sufficient  sensitivity  to  detect  a shear 
wave  in  unconsolidated  materials. 

In  the  absence  of  external  pressure,  unconsoli- 
dated sediments  display  large  shear  wave  attenuations  which 
increase  rapidly  with  frequency.  Relatively  low 
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frequencies  are  required  if  the  signal-to-noise  ratio  is  to 
remain  tolerable  over  even  short  propagation  distances. 

The  resonance  frequency  of  a ceramic  shear  element  is  high 
unless  the  length  to  thickness  ratio  is  small.  Additional- 
ly, the  ceramic  material  is  not  nearly  as  compliant  as  the 
sediment.  A large  impedance  mismatch  exists  between  the 
transducer  ,and  the  sample.  This  limits  the  amount  of  ener- 
gy transferred  to  the  sediment.  Most  of  the  laboratory 
pulse  experiments  reported  in  the  literature  for  shear  wave 
propagation  in  unconsolidated  materials  have  required  the 
application  of  external  frame  pressure  to  increase  the  ri- 
gidity to  a point  where  a signal  could  be  detected.  Al- 
though such  pressures  may  exist  at  depth  in  a sediment  col- 
umn, they  are  not  representative  of  conditions  at  the  sedi- 
ment  surface.  Comparison  of  experimental  results  therefore 
requires  specification  of  the  overburden  pressure.  Values 
of  shear  wave  speed  and  attenuation  obtained  with  minimal 
external  pressure  can  be  quite  different  from  measurements 
where  the  vertical  stress  is  due  only  to  self-weight. 

Shirley  and  Anderson  (1975a)  conducted  a series 
of  tests  designed  to  develop  a transducer  capable  of  making 
measurements  in  materials  with  shear  moduli  in  the  range  of 
natural  ocean  sediments.  Continued  development  (Shirley 


Resonant  Column  Techniques 


The  forced  vibration  of  a cylindrical  rod  is  a 
common  laboratory  means  of  determining  dynamic  properties 
at  low  frequencies.  Speeds  can  be  calculated  from  the 
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length  of  the  sample  and  the  resonance  frequency.  Both  the 
shape  of  the  resonances  in  steady  state  motion  and  the  de- 
cay of  the  amplitude  with  time  for  free  vibration  yield  a 
measure  of  the  internal  damping. 

If  the  length  of  the  cylinder  is  large  compared 
to  its  diameter,  both  torsional  and  longitudinal  vibrations 
are  approximately  one  dimensional.  For  a perfectly  elastic 
rod  fixed  at  one  end  and  free  at  the  other,  the  resonance 
frequencies  depend  on  the  length  L,  through 


u>n  = (n/2L)  /c/p  n = l,3,5,-.-» 


II-l 


where  c is  respectively  the  shear  modulus  and  Young's 
modulus  for  torsional  and  longitudinal  vibrations  (Wu  1971* 
p.  127).  The  integer  n gives  the  number  of  1/4  wavelengths 

V 

within  the  cylinder.  The  expression  will  change  for  differ- 
ent boundary  conditions. 

If  the  damping  is  small,  the  quality  factor  Q is 
obtained  from  the  steady  state  amplitude  versus  frequency 
curve  by 


Q = (w  /Au>)  ( AmRir-A2 ) /A2  » 

n max 


II-2 


where  A is  the  amplitude  at  u and  A is  the  amplitude  at 
max  1 n 

l/2w  on  either  side  of  the  resonance  (Hall  and  Richart, 
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1963).  Conventionally, 


Q = f/Af  , 


.707  A. 


max 


ment  6 . 

The  attenuation  a in  nepers/unit  distance  mea- 
sured in  pulse  experiments  is  related  to  the  decay  per 
wavelength  6 and  the  sharpness  of  the  resonance  Q,  for 
small  damping,  by 


= aV/f  = 7r / Q , 


II-4 


where  V is  the  velocity  of  the  appropriate  wave  type.  The 
attenuation  for  shear  waves  in  the  bulk  medium  and  the  de- 
cay of  torsional  vibrations  of  a cylinder  are  directly  re- 
lated by  Eq.  II-U  when  V is  the  shear  velocity.  However, 
for  longitudinal  vibrations,  6 and  Q as  determined  by  reso 
nant  column  experiments  strictly  refer  to  Young's  modulus 
waves  rather  than  the  bulk  compressional  wave  measured  by 
pulse  transmission.  This  distinction  is  not  always  made 
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when  comparing  pulse  measurements  with  resonant  column  re- 
sults. In  the  present  report,  dilational  wave  parameters 
referring  to  one  dimensional  Young's  modulus  waves  are 
termed  longitudinal,  whereas  hulk  dilational  waves  are 
called  compressional . The  differences  are  often  negligible 
(White,  1965). 

The  wave  equations  for  a perfectly  elastic  medium 
are  based  on  the  assumption  that  the  stress  can  be  ex- 
pressed as  a linear  combination  of  the  strains.  If  the  me- 
dium is  driven  by  a sinusoidal  force,  then  the  stress  and 
strain  will  be  in  phase.  In  a lossy  medium,  a steady  state 
oscillation  will  still  result  in  a stress  that  is  propor- 
tional to  the  strain,  but  the  two  will  be  out  of  phase 
(White,  1965).  The  shape  of  the  closed  curve  representing 
a stress-strain  cycle  can  thus  be  used  to  measure  the  damp- 
ing properties  of  the  material. 

In  the  case  of  a simple  shear  stress  cycle,  the 
points  of  maximum  stress  define  a line  whose  slope  is  the 
shear  modulus.  The  area  within  the  loop  is  proportional  to 
the  energy  dissipated  in  the  cycle  and  is  related  to  the 
maximum  stored  energy  by  the  damping  ratio  D, 

D = A/2*Am 


II-5 
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where  A is  the  area  of  the  loop  and  A is  the  area  beneath 

m 

the  line  defining  the  shear  modulus  (Hardin  and  Drnevich, 

1972).  For  small  damping,  6 = 2irD. 

The  condition  of  simple  shear  stress  is  obtained 
in  the  laboratory  by  application  of  a reversible  torque  to 
the  end  of  a hollow  cylinder  of  the  material  under  investi- 
gation. The  ratio  of  the  wall  thickness  to  the  radius  of 
the  cylinder  must  be  small  enough  that  radial  variations  of 
stress  are  insignificant.  The  cylindrical  shape  of  the 
specimen  avoids  problems  with  boundary  conditions  which 
would  exist  if  a shear  stress  was  imposed  on  a finite  block 
of  material  (Hardin  and  Drnevich,  1972). 

The  application  of  the  techniques  discussed  above 
present  several  practical  problems  for  unconsolidated  mate-  ^ 
rials.  Since  sediments  lack  sufficient  cohesion  to  main- 
tain themselves  in  cylindrical  shapes,  a container  must  be 
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used.  The  resonance  frequencies  will  depend  not  only  on  the 
sediment  properties  but  also  on  those  of  the  container. 

Additionally,  the  packing  of  particles  in  such  a cylinder 
may  differ  from  packings  in  pulse  experiments  on  loose  ma- 
terial, Analysis  of  the  boundary  conditions  is  involved 
and  will  vary  from  one  apparatus  to  another.  Often  the  ef- 


fects of  external  stress  are  studied.  This  causes 


additional  modification  of  the  simple  cases  discussed 
above.  The  reader  may  refer  to  the  articles  cited  or  to 
any  of  the  many  texts  on  soil  mechanics  if  additional  de- 
tail is  desired,  e.g.,  Wu  (1971)  - 

Since  investigators  may  report  the  damping  prop- 
erties of  a material  in  terms  of  a,  6,  Q,  or  D,  the  rela- 
tionships between  these  quantities  have  been  summarized  in 
Fig.  II-l.  The  formulas  are  valid  only  if  the  damping  is 
small.  Small  will  be  defined  in  a later  section. 

Laboratory  Simulation  of 
In  Situ  Conditions 

Sediments  in  the  field  are  often  subject  to 
stresses  not  found  in  small  laboratory  samples.  The  weight 
of  overlying  material  will  cause  the  rigidity  to  vary  with 
depth.  Sediments  at  the  surface  of  the  ocean  floor  are 
subjected  to  hydrostatic  pressure  not  found  in  a core  sam- 
ple. Much  of  the  work  reported  in  the  literature  has  been 
concerned  with  the  variation  of  acoustic  parameters  with 
stress  in  an  effort  to  simulate  in  situ  conditions.  These 
studies  often  give  an  acoustic  property,  such  as  the  speed, 
as  a function  of  pressure.  However,  since  stress  is  a ten- 
sor inadequately  defined  by  the  term  pressure,  the 
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different  experimental  means  of  applying  external  compac- 
tion must  be  examined  before  their  results  are  compared.  A 
precise  comparison  is  often  not  possible  without  detailed 
analysis . 

A typical  composite  of  a compression  device  can 
apply  three  different  forces  to  a cylindrical  sample.  The 
fluid  in  the  pore  space  is  connected  to  an  external  hydrau- 
lic system  that  regulates  the  pore  pressure.  A piston 
aligned  with  the  axis  of  the  cylinder  exerts  a vertical 
stress  on  the  frame  and  fluid  which  is  opposed  laterally  by 
the  confining  pressure.  The  confining  pressure  is  supplied 
by  the  container  walls  or  controlled  by  a second  hydraulic 
system  isolated  from  the  pore  space  by  a flexible  membrane. 

The  compression  of  the  sample  in  a particular  di- 
rection  is  determined  by  the  total  stress  in  that  direction 
minus  the  contribution  of  the  pore  pressure.  The  effective 
stress  is  often  defined  as  the  difference  of  the  two  (Wu, 
1971»  p.  13).  It  is  the  effective  stress  that  is  normally 
reported  as  the  variable  in  dynamic  experiments.  However, 
some  experimenters  report  the  effective  vertical  stress, 
often  called  the  differential  pressure,  while  others  re- 
port the  effective  confining  pressure.  The  two  are  not 
equivalent  (Hamilton,  1976c).  The  relationship  between  the 


two  depends  on  the  elastic  moduli  of  the  material  and  on 
the  experimental  control  of  both.  No  effort  will  be  made 
in  the  present  report  to  connect  the  two. 

Several  additional  terms  arise  in  the  discussion 
of  stress  within  natural  sediments.  For  uniform  layers  of 
sediment,  the  total  vertical  stress  at  a given  depth, 
called  the  overburden  or  geostatic  pressure,  is  determined 
by  the  weight  of  the  sediment  and  fluid  above  that  depth. 
The  hydrostatic  or  pore  pressure  is  given  by  the  weight  of 
the  water  column.  The  effective  stress  obtained  from  the 
overburden  pressure  minus  the  pore  pressure  is  sometimes 
called  the  frame  or  skeleton  pressure  and  corresponds  to 
the  differential  pressure  applied  in  the  laboratory. 

The  above  discussion  implies  that  a given  parame- 
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ter  dependent  on  the  compression  is  a constant  for  a par- 
ticular effective  pressure  regardless  of  the  combinations 
of  total  stress  and  pore  pressure  involved.  This  is 
strictly  true  only  for  special  cases.  Modifications  to  the 
definition  of  effective  stress  given  above  are  sometimes 
employed  (Robin,  1973). 

The  laboratory  control  of  physical  conditions 
such  as  the  pressure  makes  possible  a detailed  study  of  the 
effects  of  various  parameters  on  acoustic  propagation  in 


porous  media.  However,  even  when  environmental  conditions 
are  simulated  as  closely  as  possible  in  the  laboratory,  the 
exact  duplication  of  in  situ  characteristics  is  doubtful. 

In  situ  measurements  are  therefore  required  for  a complete 
acoustic  description  of  natural  unconsolidated  materials. 
Additionally,  lower  frequencies  can  be  used  in  situ  than  are 
possible  in  laboratory  pulse  techniques. 


In  Situ  Pulse  Transmission 


The  most  direct  measurements  of  in  situ  acoustic 
properties  are  obtained  from  extensions  of  the  transmission 
techniques  discussed  above.  A signal  can  be  sent  between 
probes  inserted  to  shallow  depths  in  loose  sediments.  This 
has  been  done  in  the  ocean  hy  divers  at  shallow  water 
depths  and  from  research  submersibles  approximately  1,000  m 
below  the  ocean  surface  (Hamilton  et  al.,  1970).  Addi- 
tionally, the  ocean  bottom  has  been  probed  by  transducers 
attached  to  the  cutting  head  of  a geologic  corer  (Shirley 
and  Anderson,  1975c).  The  propagation  path  is  as  near  the 
end  as  possible  to  minimize  disturbance  of  the  sediment. 


Signals  have  also  been  sent  between  two  or  more 


boreholes  , 


between  the  surface  and  a source  or  receiver  in  a single 
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borehole,  or  between  source  and  receiver  in  the  same  bore- 
hole. Since  the  propagation  path  lengths  are  often  orders 
of  magnitude  greater  than  those  in  the  laboratory,  the  com- 
pressional  wave  sources  for  large  scale  transmission  mea- 
surements are  normally  explosive  charges  or  other  high 
energy  devices.  The  detectors  are  often  moving  coil  geo- 
phones sensitive  to  velocity  changes  rather  than  pressure. 

Accurate  compressional  wave  measurements  are  eas- 
ily obtained  by  the  above  procedures,  although  special  cor- 
rections must  be  applied  for  certain  geometries.  The  gen- 
eration and  detection  of  shear  waves,  however,  is  sometimes 
more  difficult.  Shear  waves  are  normally  excited  during  an 
explosion  (Geyer  and  Martner , 1969)  and  can  be  enhanced  by 
special  shaping  of  the  charge.  Ohta  and  Shima  (1967)  found 
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a charge  detonated  in  an  iron  pipe  to  be  an  effective  gen- 
erator of  S waves.  Warrick  (197*0  tested  Primacord  on  the 
side  of  a shallow  pit.  However,  these  techniques  also  gen- 
erate compressional  waves  which,  due  to  greater  speed,  of- 
ten obscure  the  shear  wave  arrival.  If  both  source  and  re- 
ceiver are  near  the  surface,  additional  interference  will 
be  caused  by  surface  waves  whose  speeds  are  slightly  less 


than  those  of  the  shear  wave. 
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In  their  study  of  a contained  explosion,  Ohta  and 
Shima  (1976)  found  that  vertically  polarized  shear  waves 
could  easily  be  generated  by  striking  the  casing  pipe  of  a 
shot  hole  with  a hammer.  Stokoe  and  Woods  (1972)  used  a 
hammer  blow  on  a rod  inserted  in  a shallow  borehole  as  a 
shear  wave  source  for  small  scale  in  situ  tests.  White  and 
Sengbush  (1953)  found  that  a weight  dropped  into  a contain- 
er spiked  to  the  borehole  wall  was  an  adequate  generator  of 
vertically  polarized  shear  waves. 

Horizontally  polarized  shear  waves  can  be  gener- 
ated by  a blow  to  the  end  of  a wooden  beam  held  in  place  by 
the  weight  of  a truck  (Kudo  and  Shima,  1970;  Warrick, 

197*+).  A greater  impulse  can  be  applied  by  shooting  water 

or  heavy  weights  from  an  iron  pipe  fixed  to  the  beam  (Jol- 
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ly , 1956;  Shima  ^ariKd  Oiita,  1967).  If  the  direction  of  the 
blow  applied  to  the  beam  is  reversed,  the  received  shear 
wave  pulse  will  change  polarity.  However,  compressional 
wave  noise  will  not.  Thus  the  received  signal  can  be  en- 
hanced by  recording  signals  from  oppositely  directed  im- 
pulses and  subtracting  one  from  the  other.  Mechanical  vi- 
brators can  be  used  to  impart  lateral  motion  to  a plate  for 
the  production  of  shear  waves  at  a given  frequency  (Cherry 
and  Water,  1968). 
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Indirect  measurement  of  shear  properties  has 
sometimes  been  more  feasible  than  the  design  of  a shear 
wave  source.  Cunny  and  Fry  (1973)  derived  shear  wave 
speeds  from  Rayleigh  wave  velocities.  Hamilton  et  al . 
(1970)  did  the  same  for  Stoneley  waves  on  the  ocean  floor. 
The  effective  depth  of  such  an  estimate  is  approximately 
one  half  the  wave  length  of  the  boundary  wave. 

In  Situ  Reflection/Refraction 
Profiles 

In  situ  placement  of  measuring  devices  at  great 
depths,  either  in  the  sea  or  beneath  the  sediment  surface, 
is  not  economically  viable  for  surveys  of  large  areal  ex- 
tent. However,  routine  surveys  are  conducted  by  standard 
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reflection  and  refraction  techniques  that  yield  profiles 
tens  of  kilometers  in  length  and  several  kilometers  in 
depth.  The  typical  exploration  survey  is  a vertical  pro- 
file with  signals  near  normal  incidence.  Variable  incident 
angle  records  are  also  commonly  obtained  by  the  use  of  ex- 
pendable sonobuoys.  A detailed  analysis  of  these  tech- 
niques is  beyond  the  scope  of  this  paper,  but  a review  of 
the  main  features  is  necessary  for  an  evaluation  of  the 


data  in  the  literature.  The  reader  who  desires  further 


detail  may  consult  one  of  the  texts  on  exploration  geo- 
physics (e.g.,  Dobrin,  1 9T 6 ) . A summary  of  the  techniques 
as  applied  to  marine  geophysics  is  given  by  several  arti- 
cles in  The  Sea,  Vol.  3.  (Hill,  ed.  , 1963).  Also  useful  are 
papers  by  Le  Pichon  et  al.  (1968)  which  discusses  both  verti- 
cal and  variable  angle  profiling  and  by  Rutherford  (1976) 
which  emphasizes  analytical  techniques  for  determining  ve- 
locity profiles  in  unconsolidated  sediments. 

Compressional  wave  speeds  derived  from  reflection 
returns,  either  by  simple  geometric  techniques  (Grant  and 
West,  1965)  or  by  detailed  analysis  of  multifold  data  (Tan- 
er  and  Koehler,  1969),  refer  to  average  interval  velocities 
between  discontinuities.  Such  values  are  accurate  to  ap- 
proximately 100  m/sec  in  ocean  sediment  surveys  ( Le  Pichon 
et  al . , 1968).  For  deep  sea  sediments  the  earliest  returns 
are  usually  50  to  100  m below  the  sediment  surface  (Hamil- 
ton et  al . , 1977),  The  results  are  useful  for  determining 
the  structural  relationships  of  sediment  layers  and  for  es- 
tablishing velocity  profiles  with  depth.  They  do  not  yield 
detailed  information  on  the  surface  acoustic  properties. 
Unconsolidated  sediments  are  found  to  depths  over  a 1,000  m 


(Houtz  and  Ewing,  1963). 
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Although  deep  sea  sediments  do  not  return  signals 
of  significant  strength  from  the  sediment  surface,  it  is 
sometimes  possible  to  relate  an  observed  phase  change  with 
the  surface  velocity  (Fry  and  Raitt , 1961).  This  model  as- 
sumes a low  sediment  velocity  and  treats  the  sediment  as  a 
fluid. 

Conventional  exploration  reflection  surveys  pro- 
vide only  the  gross  properties  of  bottom  sediments.  How- 
ever, this  is  due  in  part  to  intentional  design  specifica- 
tions. That  is,  the  normal  seismic  profiling  system  is  de- 
signed to  outline  the  structure  of  reflective  surfaces  at 
great  depths,  not  near  the  ocean  floor.  Since  the  trans- 
mission paths  are  large  energy  losses  due  to  attenuation 
must  be  minimized.  This  requires  low  frequency,  long  wave- 
length signals  incapable  of  resolving  thin  layers.  For 
ease  of  operation,  the  source  and  receiver  are  near  the 
top  of  the  water  column.  Conversely,  the  ideal  system  for 
study  of  the  top  sediment  layers  would  include  a high  fre- 
quency source  towed  near  the  bottom.  Such  devices  exist 
(Tyce,  1976).  As  expected,  they  provide  excellent  defini- 
tion of  the  sediment  fine  structure.  It  is  also  possible 
to  obtain  estimates  of  the  speed  and  attenuation  if  certain 
simplifying  assumptions  are  made. 
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Acoustic  energy  can  be  returned  from  the  subsur- 
face by  mechanisms  other  than  simple  reflections.  The 


classical  refraction  method  of  seismic  exploration  is  based 
on  the  return  of  energy  from  a head  wave  traveling  at  graz- 
ing incidence  to  a boundary  between  isovelocity  layers. 

The  wave  speeds  and  depths  of  these  layers  can  be  obtained 
from  simple  geometric  considerations  provided  the  interval 
velocities  increase  with  increasing  depth.  In  some  appli- 
cations, the  refraction  returns  arise  from  a turning  point 
associated  with  a continuous  velocity  gradient  rather  than 
from  a boundary  wave  between  layers.  Curved  ray  path  mod- 
els based  on  assumed  velocity  gradients  allow  the  veloci- 
ties at  greatest  penetration  to  be  inferred.  Generally, 
refraction  surveys  provide  better  velocity  information  but 
less  structural  definition  than  do  reflection  profiles. 

A conventional  refraction  profile  contains  both 
direct  arrivals  and  head  wave  returns.  At  short  source-re- 
ceiver offsets  the  direct  wave  will  arrive  first.  At  larg- 
er offsets,  the  head  wave,  which  travels  at  the  greater 
speed  of  the  lower  medium,  will  overtake  the  direct  wave. 
Identification  of  the  head  wave  is  difficult  until  this 
crossover  distance  is  reached.  The  method  is  most  useful 
when  the  source  to  receiver  distance  is  large  compared  to 
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the  depth  of  the  discontinuity.  It  is  an  excellent  tech- 
nique for  obtaining  interval  velocities  to  shallow  depths 
on  land.  It  will  not,  however,  reveal  a low  velocity  lay- 
er. Since  marine  sediments  often  have  speeds  lower  than 
the  overlying  water,  conventional  refraction  techniques 
will  not  sample  the  surface  of  the  deep  ocean  floor.  For 
greater  sediment  depths,  the  refraction  record  will  be  dom- 
inated by  the  reflections  unless  the  offset  distance  is 
large  compared  to  the  depth  below  the  water  surface  or  the 
velocity  contrasts  are  small  (Kennet  and  Orcutt,  1976). 

The  interpretation  of  head  wave  returns  is  based 
on  the  assumption  of  isovelocity  plane  layers.  If  the 
planes  are  dipping,  profiles  in  two  directions  are  required 
for  correct  interpretation.  If  the  layer  velocity  is  an 
increasing  function  of  depth,  then  the  travel  path  between 
source  and  discontinuity  will  be  curved.  Analysis  in  such 
a case  is  possible  (Grant  and  Vest,  1965).  However,  it  is 
not  always  possible  to  differentiate  between  a straight 
line  head  wave  record  and  one  that  is  slightly  curved  due 
to  a moderate  velocity  gradient  (Kenneth  and  Orcutt,  1976). 

It  is  theoretically  possible  to  use  curved  path 
refractions  to  determine  the  velocity  profile  within  a sin- 
gle layer  with  a positive  velocity  gradient  (Hanna,  1973). 
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However,  refractions  near  grazing  incidence  are  hard  to 


separate  from  the  surface  reflections.  Houtz  (1978)  pre- 
sents a method  for  travel  time  inversion  of  water  multiple 
refractions  that  is  useful  in  certain  shallow  water  envi- 
ronments and  postulates  that  attenuation  can  also  be  suc- 
cessfully measured.  Case  studies  show  that  the  technique 
is  not  reliable  at  sediment  depths  of  less  than  200  m. 

Shear  wave  reflection  and  refraction  profiles  are 
similar  to  compres s ional  wave  records  (Jolly,  1956).  A 
simple  traction  source  as  discussed  above  is  adequate  for 
near-surface  records  on  land.  Multichannel  shear  wave 
seismograms  comparable  in  quality  to  compressional  wave  pro- 
files can  be  obtained  with  vibrator  arrays  as  sources.  The 
method  however  is  not  commonly  used.  It  is  also  possible 
to  observe  shear  wave  events  that  arise  from  compressional 
wave  conversion  at  an  interface.  Ricker  and  Lynn  (1950) 
delineated  subsurface  structure  from  observations  of  waves 
that  traveled  as  compressional  waves  to  the  interface  and 
returned  as  shear  waves  after  reflection. 

Since  transverse  oscillations  are  not  supported 
by  water,  shear  wave  determinations  from  marine  profiles 
must  be  made  on  the  basis  of  converted  waves.  Compression- 


al wave  arrivals  that  travel  as  shear  waves  in  the  ocean 
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floor  have  been  observed  and  verified  (Helmberger  and  Mor- 
ris, 1970).  Recent  reports  have  studied  the  interpretative 
possibilities  of  these  events  (Tatham  and  Stoffa,  1976; 
Wood,  1978).  However,  they  have  not  yet  proven  useful  for 
the  determination  of  sediment  properties. 

The  reflection  coefficient  of  marine  sediments  as 
a function  of  incident  angle  can  be  computed  from  theoreti- 
cal models.  A detailed  match  with  observations  is  possible 
only  when  finite  sediment  rigidity  is  assumed.  Although 
the  agreement  is  often  excellent,  the  inversion  of  the  pro- 
cess to  infer  sediment  properties  from  synthetic  records  is 
expensive  and  inefficient  as  well  as  non-unique  (Kennett, 
1976)  . 


Compressional  Wave  Measurements 

Although  the  present  research  is  mainly  concerned 
with  shear  wave  propagation,  it  is  necessary  10  review 
briefly  the  established  features  of  compressional  wave 
propagation  in  porous  media  in  order  to  understand  restric- 
tions placed  upon  the  various  theoretical  models  that  in- 
clude both  types  of  waves. 
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Speeds 


The  absolute  value  of  the  compressional  wave 
speed  in  a given  medium  depends  upon  the  elastic  moduli, 
particularly  on  the  compressibility.  Porous  sedimentary 
rocks  have  speeds  ranging  from  roughly  3,000  to  8,000  m/sec 
at  zero  confining  pressure  (Gregory,  1975).  Unconsoli- 
dated, water  saturated  sands  have  typical  speeds  of  approx- 
imately 1,700  m/sec.  Water  itself  has  a compressional  wave 
speed  near  1,500  m/sec.  High  porosity  saturated  clays  have 
been  found  to  have  speeds  a few  percent  lower  than  water. 
Uncompacted  dry  sands  have  speeds  as  low  as  100  m/sec. 

Once  typical  values  of  compressional  wave  speed 
have  been  established  for  the  various  media,  the  variation 
of  the  speed  with  a change ;of  one  or  more  parameters  can  be 
studied.  Parameters  investigated  in  the  literature  include 
pore  fluid,  frequency,  grain. size,  porosity,  density,  pres- 
sure, and  temperature. 

The  variation  in  speed  caused  by  a change  in  pore 
fluid  has  been  of  considerable  interest  to  the  petroleum 
industry.  Several  of  the  references  cited  in  this  paper, 
e.g.,  Wyllie  et  al . (1962),  include  measurements  with  vari- 
ous organic  liquids.  The  observed  speeds  depend  both  on 
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the  density  and  compressibility  of  the  fluids  (Brown  and 
Korringa,  1975).  In  general,  only  results  for  water,  in- 
cluding brine,  and  air  or  other  gases  will  be  reviewed 
here.  The  results  for  water  are  most  applicable  for  natur- 
al unconsolidated  marine  sediments.  Results  for  dry  mate- 
rials are  included  for  comparison  of  theoretical  predic- 
tions attendant  upon  liquid  saturation.  Three-phase  sys- 
tems, that  is,  a combination  of  gas  and  liquid  in  the 
pores,  have  also  been  extensively  studied  but  will  not  be 
considered  in  the  present  paper.  The  introduction  of  even 
small  amounts  of  air  into  an  otherwise  saturated  sediment 
can  cause  pronounced  changes  in  acoustic  propagation. 

There  is  little  evidence  for  compressional  wave 
dispersion  in  either  dry  or  saturated  porous  materials. 

V 

Schmidt  (195*0  found  a constant  longitudinal  speed  for  dry 
sand  samples  from  30  to  500  Hz.  His  resolution  appears  to 
have  been  approximately  10  percent  of  his  reported  speeds. 
Hunter  et  al . (1961)  measured  compressional  wave  speeds  in 
three  well  sorted  dry  sands.  They  note  a maximum  deviation 
of  U percent  from  the  mean  value  for  measurements  from  7 to 
73  kHz.  Neither  Nolle  et  al.  (1963)  nor  McLeroy  and 
DeLoach  (1968)  report  any  variation  of  compressional  wave 
speed  for  respective  ranges  of  U00  to  1,000  kHz  in 
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saturated  sands  and  15  to  1,500  kHz  in  natural  marine  sedi- 
ments. However,  neither  paper  presents  the  original  data 
nor  resolution  limits.  Hampton  ( 1 9 67 ) reports  a change  of 
U to  6 percent  of  the  velocity  of  water  in  fine  grained 
saturated  sediments  between  3 and  200  kHz.  Hamilton  (1972) 
was  unable  to  measure  dispersion  in  situ  over  the  narrow 
range  of  3.5  to  l4  kHz. 

Grain  size  appears  to  have  little  effect  on  the 
compressional  wave  speed.  Hunter  et  al . (1961)  noted  a 
slight  increase  in  speed  with  an  increase  in  particle  size 
for  three  samples  of  dry  sand  subject  to  a slight  frame 
pressure.  Wyllie  et  al . (1956)  found  a constant  speed  of 
1,850  m/sec  for  seven  samples  of  saturated  glass  beads 
ranging  in  size  from  .028  to  .6  mm.  Nolle  et  al . (1963) 
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observed  no  change  in  a compressional  wave  speed  of 
1,7^0  m/sec  for  four  samples  of  well  sorted  sand  with  dif- 
ferent mean  diameters. 

In  the  above  studies,  the  porosity  was  essential- 
ly constant  regardless  of  the  grain  size.  This  is  expected 
for  uniform  spherical  particles  such  as  glass  beads  since 

the  volume  of  a sphere  inscribed  in  a cube  remains  in  con- 

/ 

stant  proportion  to  the  volume  of  the  cube  regardless  of 
the  sphere  radius.  Rogers  and  Head  (l96l)  have  shown  that 
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the  porosity  of  natural  sands  is  similarly  independent  of 
grain  size  for  well  sorted  samples.  However,  differences 
in  porosity  with  grain  size  are  observed  in  natural  sedi- 
ments due  to  differences  in  compaction,  sorting,  and  grain 
shape . 

Sutton  et  al.  (1957)  present  an  empirical  formula 
for  sound  speeds  taken  from  cores  that  includes  grain  size, 
porosity,  and  percent  carbonate.  The  coefficient  of  the 
porosity  term  is  an  order  of  magnitude  greater  than  that  of 
the  grain  size.  The  relative  amount  of  carbonate  is  a mi- 
nor correction.  Brandt  (1959)  compares  the  influence  of 
various  parameters  on  compressional  wave  speed  and  con- 
cludes on  the  basis  of  data  presented  by  Hamilton  et  al . 
(1956)  that  sediment  porosity  is  the  most  important  factor 
affecting  velocities  in  water  saturated  sands.  Akal  (1972) 
plots  8,287  values  of  compressional  wave  speed  for  marine 
sediments  and  finds  a range  of  velocities  relative  to  water 
from  1.30  to  .95  for  porosities  from  .25  to  .90.  He  con- 
cludes that  not  only  is  porosity  the  most  important  factor 
influencing  the  speed  but  also  that  variations  in  the  den- 
sity of  natural  sediments  are  closely  correlated  to  the 
porosity.  Most  of  15,12U  samples  lay  between  closely 
spaced  lines  representing  linear  variations  of  density  with 


porosity  based  on  grain  densities  of  2.60  and  2.75  g/cm3. 
The  density  of  the  common  minerals  are  within  this  range 
( Hurlbut , 1971). 

The  above  reports  compare  velocity  variations 
with  porosity  among  sediment  samples  of  different  composi- 
tion. It  is  also  possible  to  study  porosity  changes  within 
a single  specimen.  This  removes  the  effects  of  grain  size 
or  shape.  Hardin  and  Richart  (1963)  studied  porosity  vari- 
ations in  several  sorted  sands  with  different  packings. 
Shirley  (1977)  reports  speed  variations  with  porosity  for 
three  different  sands,  including  two  of  the  samples  con- 
sidered in  the  present  study.  He  notes  a uniform  change  in 
speed  from  1,750  to  1,680  m/sec  for  a porosity  change  from 
.38  to  .45  in  a fine  Panama  City  beach  sand. 

V 

Velocity-density  relationships  are  often  given  as 
a means  of  determining  acoustic  impedance  (Hamilton,  1978). 
However,  since  porosity  and  density  are  linearly  related 
for  unconsolidated  marine  sediments,  the  two  parameters  re- 
present essentially  the  same  variable. 

Compaction  of  a porous  material  by  application  of 
external  pressure  can  result  in  significant  changes  of  po- 
rosity, particularly  for  clays,  which  influence  the  density 
and  wave  speed.  However,  external  stresses  also  influence 
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the  elastic  moduli  of  the  sediment  components.  This  leads 
to  velocity  variations  independent  of  the  porosity.  Sands 
subject  to  pressures  encountered  at  shallow  depths  show  lit- 
tle change  in  porosity  with  compaction  (Hamilton,  1976c). 

The  incremental  change  in  compressional  wave 
speed  for  a given  change  in  pressure  has  been  observed  to 
decrease  as  the  total  pressure  in  increased.  This  change 
is  often  represented  by  a simple  power  equation,  V = kPn . 
However,  there  is  considerable  variation  in  the  value  of 

f 

the  exponent. 

Laughton  (1957)  found  the  compressional  wave 
speed  for  fine  grained  marine  sediments  to  vary  approxi- 
mately as  the  square  root  of  the  frame  pressure  to  100  MPa 
with  atmospheric  pore  pressure  (l  MPa=106  N/m2=l45  psi). 
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He  found  only  minor  changes  of  speed  with  changes  in  fluid 
pressure  to  Uo  MPa  at  constant  frame  pressure.  Hamilton 
(1976c)  quotes  a value  of  n = .013  supplied  by  Shell  De- 
velopment Company  for  longitudinal  waves  in  fine  St.  Pet- 
er's sand  from  .lU  to  2.8  MPa.  Domenico  (1977)  found  the 
compressional  wave  speed  in  brine  saturated  Ottawa  sand  and 
glass  beads  to  be  proportional  to  the  .056  power  of  differ- 
ential pressure  to  35  MPa.  This  represents  a change  for 
the  Ottawa  sand  from  1,910  m/sec  at  2.8  MPa  to  2,220  m/sec 


at  35  MPa.  Elliott  and  Wiley  (1975)  present  similar  data 
for  a coarser  Ottawa  sand  saturated  with  distilled  water. 
They  exerted  pressures  from  7 to  63  MPa.  Their  speeds  are 
noticeably  higher  than  Domenico's  at  the  same  pressure. 
However,  their  porosities  were  lower  due  to  prior  compac- 
tion. Elliott  and  Wiley  do  not  fit  a curve  to  their  data, 
but  an  exponent  of  approximately  .06  seems  to  satisfy  the 
data.  The  discrepancy  between  the  minor  variation  of  speed 
with  pressure  in  sand  and  the  large  change  in  fine  sedi- 
ments found  by  Laughton  can  be  attributed  to  a change  in 
porosity  in  Laughton's  measurements. 

The  variation  of  compres s ional  wave  speed  with 
pressure  is  greater  in  dry  than  in  saturated  sands.  The 
difference  in  speed  between  dry  and  saturated  samples  is 


therefore  most  pronounced  at  low  pressures.  Gardner  et  al. 
(196U),  Hamilton  (1976c),  and  Domenico  (1977)  all  report 
that  the  speed  for  dry  sands  and  glass  beads  varies  as  the 
l/k  power  of  the  differential  pressure.  Schmidt  (195M  was 


able  to  measure  a change  in  longitudinal  velocity  from  105 
to  170  m/sec  presumably  caused  by  a change  in  static  pres- 
sure due  to  an  increase  in  sand  thickness  from  3 to  11  cm. 
He  also  gives  values  of  98  and  113  m/sec  respectively  for 
dry  glass  bead  samples  3 and  5 cm  in  height.  Part  of  this 
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Change  may,  however,  be  caused  by  porosity  changes  induced 
by  the  boundaries  of  the  experiment.  Both  Wylie  et  al . 
(1956)  and  Hunter  et  al . (1961)  found  the  speed  in  pulse 
experiments  to  vary  until  the  distance  between  transducer 
faces  was  larger  than  approximately  100  grain  diameters. 
Wylie  et  al . estimate  a speed  of  U00  m/sec  for  dry  glass 
beads  loaded  only  by  self-weight.  When  water  saturated, 
the  same  beads  had  a velocity  of  1,850  m/sec.  The  differ- 
ence is  80  percent  of  the  saturated  speed.  For  similar 
beads  under  pressure,  Domenico  gives  values  which  differ  by 
60  percent  at  2.8  MPa  to  30  percent  at  35  MPa. 

The  pronounced  difference  between  dry  and  satu- 
rated velocities  at  low  pressures  is  not  unique  to  uncon- 
solidated materials.  The  same  behavior  to  a lesser  degree 

V 

has  been  observed  in  low  porosity  rocks  (King,  1966;  Nur 
and  Simmons,  1969). 

In  situ  measurements  relate  velocities  to  depth 
of  overburden.  An  analysis  of  well  data  by  Faust  (1951) 
indicated  a variation  in  speed  as  the  1/6  power  of  the 
depth  and  the  duration  of  burial  for  sands  and  shale.  The 
age  can  be  replaced  by  a formation  resistivity  that  depends 
on  porosity  (Grant  and  West,  1965).  The  velocity  is  not  a 
simple  function  of  depth  because  lithif ication  depends  on 
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time.  White  and  Sengbush  (1953)  note  a smooth  increase  in 
speed  with  depth  in  surface  sands  to  the  top  of  the  water 
table  that  is  consistent  with  a variation  proportional  to 
the  1/6  power.  They  measured  a horizontal  compressional 
wave  velocity  at  the  surface  of  300  m/sec. 

An  accurate  determination  of  a similar  velocity 
profile  with  depth  in  natural  marine  sediments  is  lacking. 
Normally,  interval  velocities  from  50  to  100  m below  the 
sediment  surface  are  combined  with  surface  measurements  or 
estimates  to  project  a linear  velocity  gradient  (Nafe  and 
Drake,  1957;  Hamilton  et  al . , 1977).  The  method  is  suf- 
ficient to  distinguish  between  different  sediment  environ- 
ments, but  is  unable  to  determine  a detailed  relationship 
between  overburden  pressure  and  velocity.  Normal  gradients 
are  approximately  1 to  2 m/sec/m.  These  have  been  con- 
firmed for  the  first  few  meters  by  direct  measurements  made 
while  coring  (Shirley  and  Bell,  1978). 

Precise  comparison  of  compressional  wave  veloci- 
ties in  saturated  sediments  requires  specification  of  the 
temperature  (Brandt,  I960).  Shumway  (1958)  measured  longi- 
tudinal speeds  in  five  sediments  as  a function  of  tempera- 
ture and  found  the  variation  to  be  proportional  to  that  in 
water  alone.  Values  for  a quartz  sand  ranged  from 
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1,600  m/sec  at  0°C  to  1,800  m/sec  at  70°C.  A decrease  of 
velocity  with  an  increase  of  temperature  is  sometimes  ob- 
served in  rocks.  Timur  (1977)  measured  the  effect  of  tem- 
perature on  compressional  wave  velocity  in  two  sandstones 
and  seven  carbonates  saturated  with  brine  and  subjected  to 
substantial  pore  and  frame  pressures.  He  found  an  average 
decrease  of  1.7  percent  per  100°C  temperature  change.  For 
Boise  sandstone  the  decrease  was  from  3,350  m/sec  at  20°  to 
3,260  m/sec  at  120°.  The  decrease  was  postulated  to  arise 
from  development  of  microcracks  within  the  rock  caused  by 
differential  expansion  of  the  rock  matrix. 

Attenuation 

The  variation,  of  compressional  wave  attenuation 
with  frequency,  as  a first  approximation,  can  be  expressed 
by  a simple  power  equation  a = kfn.  For  fine  grained  sedi- 
ments and  rocks  over  a wide  frequency  range,  n - 1 . For 
well  sorted  sands,  there  is  a variation  in  reported  values 
of  n which  suggests  that  a simple  power  law  may  be  an  over- 
simplification . 

Wood  and  Weston  (l 964)  measured  attenuation  be- 
tween probes  in  tidal  mud  flats.  They  found  the  attenua- 
tion to  vary  as  f1  from  4 to  50  kHz.  Cole  (1965)  compared 
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observed  reflection  returns  from  in  situ  marine  sediments 
with  theoretical  predictions  of  bottom  loss.  He  found  that 
an  assumed  f1  dependence  of  attenuation  of  frequency  gave  a 
reasonable  fit  between  .1  and  1 Khz.  Bennet  (1967)  found 
attenuation  proportional  to  f1  between  40  and  900  Hz  by 
comparison  of  surface  and  subsurface  reflection  returns 
from  a natural  wedge  of  sediments.  Tullos  and  Reid  (1969  ) 
obtained  a linear  dependence  from  50  to  U00  Hz  from  down- 
hold  shooting  through  layers  of  sediment.  McLeroy  and 
DeLoach  (1968)  measured  attenuation  in  a number  of  natural 
sediment  samples  in  the  laboratory  and  found  n = 1 from  15 
to  1,500  kHz.  McCann  and  McCann  ( 1969 ) give  f1  for  mud  and 
fl.26  for  sand  between  5 and  50  Hz.  Nolle  et  al . (1963) 
found  an  f1/2  dependence  between  U00  and  1,000  kHz  in  four 
well  sorted  sands.  Hampton  (1967)  gives  f1*37  for  clay  and 
f1/2  for  sand  from  U .to  600  kHz.  Busby  and  Richardson 
(1957)  found  the  variation  of  attenuation  with  frequency  to 
change  with  grain  size  for  several  samples  of  sand  and 
glass  beads  from  .5  to  3 MHz.  The  beads  exhibited  an  f4 
dependence.  For  these  frequencies,  the  wavelength  becomes 
comparable  to  the  grain  size.  Kuster  and  Toksoz  (197*0 
found  that  a linear  dependence  of  attenuation  on  frequency 
was  inadequate  for  measurements  in  suspensions  of  glass  and 
plastic  beads  in  various  liquids  from  100  to  500  kHz. 


k 


39 

The  values  given  above  pertain  to  measurements  in 
specific  samples  over  a fairly  wide  frequency  range.  Num- 
erous measurements  of  attenuation  at  discrete  frequencies 
in  different  sediment  samples  are  available  in  the  litera- 
ture, e.g.,  Hamilton  et  al . (1956),  Shumway  (i960).  Hamil- 
ton (1972)  has  compiled  the  data  available  from  numerous 
sources  and  presented  the  result  as  a graph  of  attenuation 
versus  frequency  for  unconsolidated  sediments  from  approxi- 
mately 1 Hz  to  1 MHz.  Even  though  the  measurements  repre- 
sent sediments  with  a wide  range  of  porosity  and  grain 
sizes,  a definite  trend  in  the  data  exists.  This  trend  is 
adequately  represented  as  an  f1  dependence  of  attenuation 
on  frequency.  However,  measurements  in  a single  sediment, 
such  as  those  above,  can  show  considerable  deviation  from 

V 

f1  behavior  over  approximately  a decade  of  frequency  and 
still  fall  within  the  scatter  of  the  overall  data.  The  de- 
tailed behavior  of  attenuation  with  frequency  over  several 
decades  is  still  not  established,  particularly  for  sand 
size  particles.  A recent  version  of  Hamilton's  graph  (Ham- 
ilton, 1976a)  is  reproduced  in  Fig.  II-2. 

Wyllie  et  al . (1962)  summarize  the  attenuation 
data  then  available  for  consolidated  sediments.  Included 
are  measurements  by  Born  (19^1)  from  1 to  4 kHz  which  imply 
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a constant  log  decrement  for  extenslonal  vibrations,  and 
measurements  by  McDonal  et  al . (1958)  in  Pierre  shale  that 
yield  an  f1  dependence  from  1 to  6 00  Hz.  Attvell  and  Ra- 
mana  (1966)  present  a plot  similar  to  Fig.  II-2  for  mea- 
surements in  rock  samples.  There  is  a single  trend  from  1 
to  107  Hz  in  agreement  with  a linear  dependence  of  attenua- 
tion with  frequency.  Considerable  scatter  exists,  however, 
when  different  rock  samples  are  compared.  Recently,  Toksoz 
et  al . (1978)  report  a linear  increase  of  attenuation  with 
frequency  from  .1  to  1.5  MHz  for  both  dry  and  saturated 
rock  samples. 

Studies  of  compressional  wave  attenuation  versus 
frequency  in  dry  unconsolidated  materials  are  not  preva- 
lent. Kottonski  and  Malechi  (1958)  measured  attenuation  to 
vary  as  f’23  from  300  dB/m  at  30  kHz  to  1,650  dB/m  at 
90  kHz.  Nyborg  et  al.  (1950)  give  values  in  sand  at  10, 

18,  and  26  kHz  which  are  not  adequately  represented  by  a 
simple  power  equation.  Pilbeam  and  Vaisnys  (1973)  found 
constant  Q values  from  1 to  20  kHz  for  glass  bead  samples 
at  low  confining  pressure. 

The  expected  magnitude  of  compressional  wave  at- 
tenuation in  surficial  marine  sediments  can  be  obtained 
from  the  data  shown  in  Fig.  II-2.  Hamilton  (1976b) 


obtained  an  average  value  of  the  log  decrement,  6 = aV/f 
for  sand  of  .10  ± .03  and  for  silt  clays  of  .02  ± .01. 
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Intuitively,  the  attenuation  in  a pure  solid  or 
pure  liquid  is  less  than  the  attenuation  in  an  aggregate 
mixture.  Shumvay  (i960)  noticed  that  attenuation  in  marine 
sediments  at  first  increases  with  increasing  porosity, 
peaks  at  roughly  0 = .5,  and  then  decreases  as  the  porosity 
continues  to  increase.  This  same  behavior  was  noted  by 
Hamilton  (1972)  from  an  analysis  of  the  data  available  in 
the  literature  for  which  porosity  values  are  given.  A 
similar  relationship  results  when  the  porosity  is  replaced 
by  the  grain  size.  In  neither  case  is  the  variation  well 
established.  Considerable  scatter  exists  when  the  overall 
data  are  compared.  An  ideal  comparison  would  separate  the 
effects  due  to  grain  size  from  those  due  to  porosity.  Evi- 
dence indicates  that  the  two  are  independent  variables. 
Nolle  et  al.  (1963)  found  the  attenuation  to  vary  as  the 
square  root  of  the  mean  grain  size  for  four  well  sorted 
sands  having  essentially  the  same  porosity.  Compressional 
wave  attenuation  has  also  been  reported  to  decrease  with 
gravitational  compaction  in  laboratory  sediments  (Shirley 
and  Bell,  1978),  that  is,  for  a constant  grain  size  and  de- 


creasing porosity. 
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Compared  with  the  study  of  pressure  effects  on 
compressional  wave  speed,  there  is  little  data  available  on 
attenuation  versus  pressure  in  unconsolidated  materials. 
Hunter  et  al.  (1961)  found  the  attenuation  in  dry  sand  to 
decrease  with  an  increase  of  frame  pressure.  Hardin  and 
Richard  (1963)  found  the  decrement  for  longitudinal  vibra- 
tions to  decrease  with  increasing  confining  pressure  for 
both  dry  and  saturated  samples  of  Ottawa  sand  and  glass 
beads.  The  damping  was  greater  in  the  saturated  samples 
than  in  the  dry.  The  measured  decrement  for  the  dry  mate- 
rials snowed  considerable  variation  with  both  pressure  and 
amplitude  of  vibration.  Gardner  et  al . (1964)  found  the 
log  decrement  in  dry  sands  to  decrease  as  the  1/6  power  of 
the  overburden  pressure.  Pilbeam  and  Vaisnys  (1973)  found 
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the  damping  in  samples  of  dry  glass  beads  to  decrease  with 
increasing  confining  pressure  from  approximately  zero  to 
one  atmosphere  as  the  .3  to  .8  power  of  the  pressure. 

Toksoz  et  al.  (1978)  and  Johnson  and  Toksoz 
(1978)  have  studied  compressional  wave  attenuation  versus 
pressure  in  porous  rocks.  They  found  the  attenuation  to 
decrease  with  increasing  differential  pressure  for  both  dry 
and  saturated  samples.  The  attenuation  was  greater  in 


saturated  than  in  dry  specimens 


The  available  data  on  shear  wave  speeds  in  fully 
saturated  unconsolidated  sediments  at  small  confining  pres- 
sure are  sparse,  both  for  in  situ  and  laboratory  measure- 
ments. Warrick  (197*0  measured  a speed  of  90  m/sec  near 
the  top  of  a bore  hole  in  a mud  flat.  Shirley  and  Anderson 
(1975)  observed  a smooth  increase  in  speeds  from  8 to 
2k  m/sec  in  kaolinite  clay  which  was  allowed  to  settle  un- 
der the  influence  of  gravity  and  interparticle  forces. 
Shirley  (1977)  reports  speeds  of  60  to  110  m/sec,  depending 
on  porosity,  for  three  sand  samples.  Shirley  and  Bell 
(1978)  measured  speeds  of  100,  U0,  and  16  m/sec,  respec- 
tively, for  sand,  silt,  and  clay  with  probes  9 cm  below  the 
sediment  surface.  The  speeds  were  observed  to  vary  as  much 
as  50  percent  with  porosity  changes  affected  by  stirring 
the  sediments. 

Davies  ( 1 9 6 5 ) calculated  shear  wave  speeds  from 
Stoneley  waves  excited  on  the  ocean  floor  by  explosive 
charges.  The  speeds  varied  from  50  to  190  m/sec  with  depth 
to  16  m.  Hamilton  et  al . (1970)  also  obtained  in  situ 
shear  wave  speeds  from  Stoneley  wave  observations.  They 
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confining  pressure  vas  exerted  by  the  atmosphere.  Hamilton 
(1976c)  reviewed  torsional  vibration  data  obtained  from 
Shell  Development  Company  for  dry  and  brine  saturated 
St.  Peters  sand.  The  shear  wave  speed  was  observed  to  vary 
as  the  1/1*  power  of  the  differential  pressure  from  approxi- 
mately 280  m/sec  at  70  kPa  to  900  m/sec  at  7 MPa.  Domenico 
(1977)  obtained  a similar  1/1*  power  variation  of  shear  wave 
speed  with  differential  pressures  from  3 to  35  MPa  for  both 
Ottawa  sand  and  glass  beads, 

Hamilton  (1976c)  recommends  a depth  exponent  of 
1/1*  for  prediction  of  shear  wave  velocities  in  sand  based 
on  available  in  situ  data.  There  is,  however,  considerable 
scatter  in  the  data  with  bounds  of  1/3  to  1/6  suggested  by 
Hamilton.  Much  of  the  data  reviewed  by  Hamilton  was  ob- 
tained  from  land  measurements.  In  some  cases  the  water  ta- 
ble was  less  than  2 m below  the  surface;  in  others  it  was 
beyond  the  depth  of  the  measurements.  Since  the  measure- 
ments were  obtained  in  natural  soils  in  the  field,  it  must 
be  assumed  that  the  sediments  above  the  water  table  ex- 
hibited various  degrees  of  partial  saturation. 

The  review  of  compressional  wave  propagation  in 
porous  media  showed  a significant  difference  between  com- 
pressional wave  speeds  in  dry  and  saturated  materials. 
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Three-phase  systems  were  not  considered,  but  data  in  the 
literature  indicate  that  major  changes  in  propagation  oc- 
cur when  only  small  amounts  of  either  gas  or  liquid  are 
added  to  an  otherwise  single  phase  pore  fluid  (e.g.,  Gard- 
ner et  al . , 196U).  In  particular,  the  variation  of  atten- 
uation with  frequency  is  definitely  nonlinear  (Anderson, 
197*0  . In  comparison,  shear  wave  speeds  show  smaller  vari- 
ations with  percent  saturation.  Hardin  and  Richart  (1963) 
found  that  the  shear  modulus  for  mature  sands  decreases  as 
much  as  15  percent  with  increasing  saturation,  with  most  of 
the  change  occurring  at  the  onset  of  liquid  introduction. 
Angular  sands  showed  smaller  variations.  Domenico  (1977) 
reports  a change  in  speed  from  691  m/sec  at  total  brine 
saturation  to  61*0  m/sec  at  total  gas  saturation  for  a dif- 
ferential pressure  of  7 MPa.  The  decrease  is  smooth  for 
intermediate  saturations.  The  error  introduced  in  Hamil- 
ton's analysis  by  consideration  of  partially  saturated 
sands  is  probably  smaller  than  the  scatter  in  the  data. 
However,  care  must  be  exercised  in  applying  the  data  to 
other  situations. 

In  each  of  the  above  studies,  i.e.,  Hardin  and 
Richart  (1963),  Hamilton  (1976c),  and  Domenico  (1977)  » the 
shear  wave  speed  was  always  less  in  the  saturated  than  in 


the  dry  material.  The  same  general  behavior  is  found  in 
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consolidated  samples  although  at  high  differential  pressure 
the  closure  of  microcracks  in  rocks  sometimes  leads  to  a 
greater  speed  in  the  dry  material  (Wyllie  et  al . , 1962; 
Gregory,  1976;  Johnson  and  Toksoz  , 1978). 

The  effects  of  grain  size  and  shape  on  shear  wave 
speeds  were  studied  by  Hardin  and  Richart  ( 1 9 6 3 ) • They 
found  the  velocity  to  be  independent  of  grain  size  and  gra- 
dation for  sands  except  for  associated  changes  in  porosity. 
Results  for  fine  grained  silts  were  dependent  on  time  and 
precompression.  Angular  sands  had  higher  speeds  than 
rounded  sands.  This  is  also  supported  by  Pilbeam  and 
Vaisnys  ( 1973 ) . 

The  variation  of  shear  wave  speed  with  frequency 

V 

in  unconsolidated  sediments  has  received  little  atten- 
tion. Schmidt  (1954)  found  no  evidence  of  dispersion  in 
dry  samples  of  sand  and  glass  beads  between  50  and  200  Hz. 
His  resolution  appears  to  have  been  within  10  percent. 

Kudo  and  Shima  (1970)  measured  a constant  phase  velocity 
between  30  and  90  Hz  in  soil.  Pilbeam  and  Vaisinys  (1973) 
found  no  variation  in  speeds  derived  from  torsional  vibra- 
tions of  dry  glass  beads  between  1-20  kHz. 


Shear  wave  velocity  measurements  are  available 
over  a larger  frequency  range  in  rocks.  Simmons  (196U) 
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compared  speeds  obtained  from  pulse  experiments  at  1-5  MHz 
with  torsional  velocities  measured  by  Birch  and  Bancroft 
(1938)  in  the  same  rock  samples  at  a few  kHz.  The  results 
differed  at  most  by  2 percent.  Peselnick  and  Outerbridge 
(1961)  used  three  different  experimental  techniques  to  mea- 
sure shear  wave  velocities  in  a sample  of  Solenhofen  lime- 
stone. No  dispersion  was  found  between  3 Hz  and  10  MHz. 

Attenuation 

Damping  in  unconsolidated  materials  has  been 
found  to  depend  on  the  confining  pressure,  the  amplitude  of 
vibration,  the  frequency,  the  degree  of  saturation,  and  the 
porosity.  There  are,  however,  much  less  data  available  on 
shear  wave  damping  than  for  compressional  waves. 

Schmidt  (195^)  found  the  shear  wave  decrement  to 
decrease  with  an  increase  of  sample  length  for  dry  sands. 
Hardin  and  Richart  (1963)  found  the  decrement  of  torsional 


vibrations  to  decrease  with  an  increase  of  confining  pres- 
sure for  both  dry  and  saturated  samples.  However,  they 
also  found  the  decrement  to  decrease  as  the  l/U  power  of 
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decreasing  vibration  amplitude.  The  amplitude  curves  were 
not  parallel  for  different  confining  pressures.  Gardner 
et  al.  (1964)  found  the  torsional  decrements  for  dry  sand 
and  glass  beads  to  decrease  as  the  1/6  power  of  the  pres- 
sure below  7 MPa  and  to  approach  a constant  value  above 
35  MPa.  Their  measurements  were  independent  of  amplitude 
at  the  levels  used.  Pilbeam  and  Vaisnys  (1973)  found  the 
torsional  decrements  for  several  sizes  of  dry  glass  beads 
to  vary  between  the  .3  and  .8  power  of  confining  pressure 
below  .1  MPa.  Their  results  were  independent  of  vibration 
amplitudes  differing  by  a power  of  10. 

Schmidt  (1954)  found  that  the  damping  was  inde- 
pendent of  frequency  for  dry  sands  between  100  and 
2,000  Hz.  Barkan  (1962)  reports  a constant  damping  ratio 
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for  soils  between  .6  and  1.8  kHz.  Kudo  and  Shima  (1970) 
made  borehole  measurements  of  attenuation  with  a surface 
shear  wave  source.  They  report  that  the  attenuation  is  ap- 
proximately proportional  to  the  first  power  of  the  frequen- 
cy from  10  to  90  Hz.  Their  results  are  reproduced  in 
Fig.  II-3.  A linear  variation  with  frequency  is  a better 
fit  for  fine  grained  soils  than  for  the  sand.  The  materi- 
als should  be  considered  as  partially  saturated.  Also, 
since  the  measurements  were  made  between  the  surface  and  a 


FREQUENCY  - H> 

DILUVIAL  SAND:  V-243»A*« 
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TERTIARY  MUSTONE:  V-420»/**« 
£ -0.5 


FREQUENCY  - H* 

ALLUVIAL  SILT:  V-102  »/••* 
£ -0.14 


KWANTOLOAM:  V-150  m/t.c 
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FIGURE  1 1-3 

INSITU  MEASUREMENTS  OF  SHEAR  WAVE  ATTENUATION  vtrius  FREQUENCY 

(after  Kudo  and  Shima  1970) 
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depth  of  UO  m,  the  results  are  averages  over  any  changes 

I 

caused  by  overburden  pressure.  Hardin  and  Drnevich  (1972) 
found  a slight  increase  in  damping  ratio  between  measure- 
ments at  .1  Hz  and  measurements  at  25  and  38  Hz  in  dry 
silt.  There  does  not  appear  to  be  any  study  of  shear  wave 


attenuation  versus  frequency  for  completely  saturated  un- 
consolidated materials.  This  is  a significant  omission 
since  damping  does  depend  on  the  degree  of  saturation. 

Hardin  and  Richart  (1963)  found  that  water  sat- 
uration increased  the  log  decrement  for  Ottawa  sand  by  a 
factor  of  1.5  to  U times  that  for  the  dry  material,  depend- 
ing on  the  amplitude  of  vibration  and  the  confining  pres- 
sure. Gardner  et  al.  (I96U)  found  the  addition  of  pore  wa- 
ter to  increase  the  decrement  from  approximately  .005  to 


.06  for  sand.  The  increase  was  not  linear  with  partial 
saturation  but  depended  on  the  sample. 

There  is  a relatively  small  amount  of  scattered 
data  from  which  to  draw  conclusions  concerning  the  effects 
of  grain  size  and  porosity  on  shear  wave  damping.  Hamilton 
(1976b)  reviews  the  available  literature  and  concludes  that 
most  values  of  the  log  decrement  corresponding  to  near-sur- 
face pressures  are  between  .1  and  .6  for  both  sa..<is  and 
clays.  He  recommends  approximate  values  for  saturated  sand 


of  .30  ± .15  and  for  saturated  silt-clays  of  .2  ± .1.  That 
is,  decrements  in  high  porosity,  fine  grained  materials  ap- 
pear to  "be  slightly  lower  than  in  relatively  low  porosity, 
coarser  sands.  Gardner  et  al . (196U)  found  the  damping  to 
he  independent  of  grain  size  and  porosity  for  sand.  Pil- 
beam  and  Vaisnys  (1973)  note  a slight  change  in  log  decre- 
ment for  different  size  dry  glass  beads.  Hardin  and 
Drnevich  (1972)  found  a decrease  in  damping  ratio  with  an 
increase  in  porosity  for  undisturbed  cohesive  soils. 

Shear  wave  attenuation  in  unconsolidated  sedi- 
ments may  be  compared  with  measurements  obtained  in  con- 
solidated materials.  In  situ  shear  wave  attenuation  mea- 
surements in  Pierre  shale  (McDonald  et  al.,  1958)  were  ap- 
proximately linear  with  frequency  from  25  to  150  Hz  al- 
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though  a definite  change  in  slope  was  apparent  from  100  to 
150  Hz.  Gardner  et  al.  (1967)  obtained  measurements  of  the 
damping  in  dry  Solenhofen  limestone  from  10  Hz  to  20  kHz. 
They  found  a constant  decrement  of  .00U3  for  all  frequen- 
cies. They  compared  this  with  a value  of  .005  at  3.59  Hz 
(Peselnick  and  Outerbridge,  1961)  and  of  .0167  at  1 MHz 
(Peselnick  and  Zietz,  1959)  obtained  in  similar  samples. 
Toksoz  et  al.  (1978)  and  Johnson  and  Toksoz  (1978)  report 


that  shear  wave  attenuation  is  linearly  proportional  to 


frequency  from  approximately  100  kHz  to  1 MHz  in  both  dry 
and  saturated  samples  of  sandstone  and  limestone.  The  at- 
tenuation decreased  with  increasing  differential  pressure 
and  was  smaller  when  the  rocks  were  dry  rather  than  sat- 
urated. 


Summary  of  Known  Properties 

The  following  salient  features  of  acoustic  prop- 
agation in  porous  media  may  be  stated. 

1.  Dispersion  of  either  compressional  or  shear 
waves,  if  present,  is  small. 

2.  Compressional  wave  speeds  in  saturated  mate- 
rials are  substantially  higher  than  in  dry  materials. 

Shear  wave  speeds  in  saturated  materials  are  slightly  lower 
than  in  dry  materials. 

3.  Compressional  wave  attenuation  varies  roughly 
as  the  first  power  of  the  frequency  over  a wide  frequency 
range,  although  for  sands  the  variation  may  differ  somewhat 
over  narrow  frequency  bands.  Shear  wave  attenuation  ap- 
pears to  follow  the  same  pattern.  There  is,  however,  a lack 
of  shear  wave  data,  particularly  for  unconsolidated,  sat- 


urated media. 


Both  compressional  and  shear  wave  speeds  in 


crease  with  increasing  effective  pressure,  with  a corres 


ponding  decrease  in  attenuation.  Shear  wave  parameters  ap 


pressures  than  do  the  compressional  wave  properties 


General  Review 

Numerous  attempts  have  been  made  to  fit  the  above 
observations  into  the  framework  of  a mathematical  theory. 
These  attempts  can  be  divided  into  two  broad  categories  de- 
pending upon  the  extent  a porous  medium  is  considered  to  be 
homogeneous.  One  approach  models  the  medium  as  a continuum 
with  viscoelastic  properties  pertaining  only  to  the  bulk 
material.  The  system  can  be  described  by  complex  moduli 
and  creep  or  relaxation  functions  chosen  to  fit  empirical 
relationships.  A second  approach  assumes  a definite  loss 
mechanism  based  on  the  interactions  of  the  constituents  of 
the  material.  This  allows  the  effects  of  observable  macro- 
scopic properties  of  the  system,  such  as  porosity  or  fluid 
compressibility,  to  be  predicted  by  the  model. 

Perhaps  the  simplest  viscoelastic  model  of  a con- 
tinuum is  obtained  by  modeling  the  substance  as  the  paral- 
lel addition  of  an  ideal  spring  and  an  ideal  dashpot . The 
combination  is  known  as  a Voigt  element  (Bland,  i960). 
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Mathematically,  this  adds  a term  proportional  to  the  rate 
of  strain  to  the  elastic,  or  Hookean,  stress-strain  rela- 
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tion.  The  resulting  equation  of  motion  is  linear  and  con- 
tains two  constants,  the  normal  elastic  modulus  and  a coef- 
ficient of  viscosity.  A harmonic  solution  at  low  frequen- 
cies yields  attenuation  proportional  to  the  square  of  the 
frequency  with  a velocity  independent  of  frequency.  At 
high  frequencies  both  the  attenuation  and  velocity  become 
proportional  to  the  square  root  of  the  frequency  (White, 
1965;  Grant  and  West,  1965), 

The  response  of  a Voigt  element  to  a constant 
stress  is  called  a creep  function  since,  in  general,  the  ex- 
tension, or  strain,  is  a function  of  time.  The  ratio  of 
the  coefficient  of  viscosity  and  the  elastic  modulus  re- 
presents  the  time  required  for  the  creep  function  to  reach 
a certain  proportion  of  its  final  value.  This  ratio  is 
known  as  the  retardation  time  (Bland,  i960).  Thus  the  co- 
efficient of  viscosity  may  be  replaced  by  a constant  with 
units  of  time.  The  simple  Voigt  model  contains  one  such 
time  constant  in  addition  to  the  elastic  modulus. 

If  the  Voigt  model  is  modified  by  an  additional 
term  proportional  to  the  rate  of  stress,  then  two  time  con- 
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stants  are  required.  Proper  selection  of  these  constants. 


allowed  Horton  (1959)  to  successfully  apply  such  a model  to 
the  attenuation  data  of  the  Pierre  shale  supplied  by 
McDonal  et  al . (1958).  Four  time  constants  and  two  moduli 
were  required  to  fit  both  the  compressional  and  shear  wave 
data  over  the  given  frequency  range.  However,  an  approxi- 
mately linear  variation  of  attenuation  with  frequency  over 
an  extended  frequency  range  would  have  required  the  addi- 
tion of  higher  order  time  derivatives  to  the  stress  strain 
equation.  These  in  turn  would  have  introduced  a predicted 
dispersion  not  in  accord  with  the  general  observations 
(White,  1965). 

If  the  number  of  Voigt  elements  added  to  the 
viscoelastic  model  increases  indefinitely,  then  it  is  pos- 
sible to  define  a distribution  function  of  retardation 
times.  Chae  (1968)  assumed  a distribution  function  for  the 
viscoelastic  properties  of  sand  which  gives  a reasonable 
fit  to  the  shear  wave  data  of  Kudo  and  Shima  (1970).  The 
latter  authors,  however,  were  able  to  obtain  slightly  bet- 
ter agreement  by  assuming  a different  distribution  func- 
tion. Neither  predicts  a strictly  constant  log  decrement 
or  phase  velocity.  Indeed,  it  has  been  shown  (Collins  and 
Lee,  1956;  Knopoff  and  MacDonald,  1958)  that  a constant 
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imaginary  and  real  parts  of  a complex  shear  modulus.  That 
is,  the  damping  is  proportional  to  tan  <|>  = y ' /y  where 
y = y + iy'.  Thus  a viscoelastic  model  results  if  the 
shear  modulus  is  regarded  as  complex.  If  the  damping  is 
small,  |y|  .«  y and  the  frequency  independent  elastic  expres- 
sion for  the  velocity,  v = (y/p)1/2,  is  recovered. 

The  above  model  is  favored  by  Hamilton  (1976b) 
since  it  is  able  to  satisfy  the  two  main  tenets  of  acoustic 
propagation  in  porous  media,  namely  that  dispersion  is  neg- 
ligible and  that  attenuation  varies  linearly  with  frequen- 
cy. Further,  the  model  indicates  how  well  the  viscoelastic 
parameters  can  be  approximated  by  the  elastic  equations. 
However,  there  is  no  provision  for  the  determination  of  the 
real  and  imaginary  parts  of  the  complex  modulus  apart  from 
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direct  measurement  of  the  speed  and  attenuation  in  the  bulk 
medium.  A theory  would  display  greater  utility  if  these 
parameters  could  be  estimated  from  the  physical  properties 
of  the  medium. 

Three  main  mechanisms  have  been  considered  as 
probable  causes  of  attenuation  in  porous  materials.  These 
are  scattering,  frictional  forces  at  grain  contacts,  and 
viscous  flow  of  pore  fluid  relative  to  the  frame.  The  to- 
tal observed  attenuation  is  probably  a combined  effect  of 


all  three  mechanisms.  However,  under  certain  conditions, 
one  or  two  mechanisms  may  dominate. 

In  a suspension,  grain  contacts  are  absent. 

Urick  (19^8)  presents  a formula  for  the  attenuation  in  a 
suspension  that  includes  both  scattering  from  the  individ- 
ual particles  and  viscous  loss  due  to  relative  motion  of 
the  particles  and  the  liquid.  At  low  frequencies,  the  vis- 
cous term  predominates;  at  high  frequencies,  scattering 
makes  the  major  contribution  (Anderson,  197M  . In  a gran- 
ular material,  scattering  is  again  important  at  high  fre- 
quencies, that  is,  when  the  wavelength  approaches  the  parti 
cle  size.  Plona  and  Tsang  (1979)  use  scattering  theory  to 
predict  that  attenuation  in  a porous  material  is  propor- 
tional to  the  cube  of  the  grain  size  and  the  fourth  power 
of  the  frequency.  They  present  experimental  confirmation 
in  saturated  glass  beads  at  frequencies  near  1 MHz.  This 
same  dependence  was  observed  by  Busby  and  Richardson 
(1957)  . 

A necessary  concept  in  the  comparison  of  contact 
losses  and  viscous  dissipation  is  the  existence  of  a 
"frame."  Wood  (1930)  developed  a formula  for  compress ional 
wave  speed  which  related  the  individual  compressiblities 
and  densities  of  the  constituents  of  a suspension  to  a bulk 
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compressibility  and  density.  The  relative  contributions  of 
the  fluid  and  the  particles  depended  on  the  porosity.  The 
Wood  equation  correctly  predicts  that  the  speed  in  high  po- 
rosity marine  sediments  will  be  lower  than  the  speed  in  wa- 
ter (Anderson,  197^).  However,  predicted  values  are  lower 
than  observed  values  (Shirley  and  Bell,  1978).  This  can  be 
attributed  to  a finite  rigidity  in  the  sediment  that  arises 
from  the  grain  to  grain  contacts.  This  is  not  the  rigidity 
of  the  material  comprising  the  grains,  which  is  quite  high, 
nor  the  rigidity  of  the  fluid,  which  has  none.  Rather,  it 
is  a property  of  the  arrangement  of  the  particles,  that  is, 
the  frame . 

Gassman  (1951)  derived  expressions  for  the  bulk 
moduli  of  a porous  material  in  terms  of  the  moduli  of  the 
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fluid,  the  grains,  and  the  frame.  The  formula  reduces  to 
that  of  Wood  (1930)  when  the  frame  moduli  are  zero.  Gass- 
man suggested  that  the  frame  moduli  could  be  measured  in 
the  absence  of  the  pore  fluid.  The  theoretical  character- 
ization of  the  frame  is  difficult. 

Gassman  studied  the  contact  areas  of  an  hexagonal 
array  of  perfect  spheres  under  the  influence  of  self- 
weight. Tangential  stress  was  neglected.  He  predicted 
that  the  speed  should  vary  as  the  1/6  power  of  the  depth. 
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Other  authors  have  considered  different  packings  of  ideal- 
ized spheres.  White  (1965)  reviews  much  of  this  work. 
Walton  (1975 » 1977)  calculates  the  moduli  of  a fluid  satu- 
rated cubic  array.  He  predicts  that  the  velocity  in  the 
dry  material  should  vary  as  the  1/3  power  of  the  depth  at 
shallow  depths  and  as  the  1/6  power  at  larger  depths.  He 
also  predicts  the  existence  of  two  compressional  waves,  one 
associated  with  the  frame  and  one  with  the  fluid. 

Attenuation  in  the  frame  can  be  introduced 
through  conjectural  frictional  forces  at  the  grain  con- 
tacts. White  (1965)  reviews  the  work,  both  theoretical  and 
experimental,  which  has  been  done  with  arrays  of  spherical 
particles.  He  concludes  that  both  static  and  dynamic  fric- 
tion must  be  considered  in  the  displacement  cycle  if  the 
predicted  attenuation  is  to  agree  with  experimental  evi- 
dence. In  consolidated  materials,  the  possibility  exists 
of  friction  between  crack  surfaces  as  well  as  at  grain 
boundaries  (Walsh,  1966).  Johnston  et  al.  (1978)  review 
this  subject.  Although  the  important  features  of  the  frame 
attenuation  can  be  predicted  qualitatively,  quantitative 
results  are  hampered  by  several  unknown  parameters.  There- 
fore the  relative  importance  of  frame  attenuation  to  vis- 
cous attenuation  must  be  based  on  observation. 
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Numerical  values  of  the  attenuation  due  to  rela- 
tive movement  of  the  fluid  and  frame  can  be  calculated. 

Biot  (1956a,  1956b)  developed  a general  theory  of  acoustic 
propagation  in  porous  media  based  on  the  earlier  vork  of 
Zwikker  and  Kosten  (1949)  and  Morse  (1952).  When  no  rela- 
tive movement  between  frame  and  fluid  exists,  Biot's  equa- 
tions reduce  to  Gassman's  (1951)  formulation  (Geertsma  and 
Smith,  1961).  When  relative  movement  is  allowed,  the  pre- 
dictions depend  upon  frequency.  At  low  frequencies,  the 
flow  in  the  pores  is  laminar,  the  speed  is  essentially  con- 
stant, and  the  attenuation  varies  as  the  square  of  the  fre- 
quency. At  high  frequencies,  the  flow  pattern  is  complex, 
the  speed  is  again  approximately  constant  although  differ- 
ent, and  the  attenuation  varies  as  the  square  root  of  the 
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frequency.  At  intermediate  frequencies,  a transition  zone 
exists.  The  details  of  the  transition  depend  upon  the  com- 
pressibility and  viscosity  of  the  fluid  and  the  moduli  and 
permeability  of  the  frame.  The  theory  assumes  a narrow 
distribution  of  pore  sizes. 

The  above  observations  apply  to  both  the  compres- 


sional  wave  and  the  shear  wave.  However,  the  theory  also 
predicts  the  existence  of  a second  dilational  wave  of  low 
velocity  and  extremely  high  attenuation.  This  second  type 
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of  compressional  wave  has  not  been  observed  in  liquid  sat- 
urated porous  media  (Yew  and  Jogi,  1976).  Even  in  materi- 
als where  the  viscous  effect  is  expected  to  be  large,  the 
attenuation  of  the  second  wave  is  expected  to  make  direct 
observation  extremely  difficult.  It  is  conceivable,  how- 
ever, that  the  presence  of  the  second  wave  could  be  deduced 
from  a change  in  reflection  coefficient  with  frequency 
since  mode  conversion  should  exist  at  a boundary.  The  sec- 
ond wave  has  been  considered  a prominent  mechanism  of  ener- 
gy loss  in  gas  saturated  acoustic  absorbers  (Zwikker  and 
Kosten,  19^9*  Ferrero  and  Sacerdote,  1951). 

The  viscous  interaction  model  correctly  predicts 
that  the  compressional  wave  speed  increases  and  the  shear 
wave  speed  decreases  with  fluid  saturation.  However,  the 
lack  of  experimental  observations  supporting  a general  f2 
or  f1/2  dependence  of  attenuation  on  frequency  and  the  ap- 
parent lack  of  dispersion  have  been  cited  as  evidence  that 
viscous  loss  is  of  minor  importance  to  overall  sediment  at- 
tenuation (Hamilton,  1972).  Although  such  a conclusion  may 
be  correct,  the  available  experimental  data  are  inadequate 
to  dismiss  the  theory  in  all  cases. 

Several  factors  must  be  considered  when  comparing 
the  theoretical  predictions  with  observations.  First  the 
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attenuation  will  be  a combination  of  both  the  contact 
losses  and  the  viscous  term.  The  contact  losses  are  diffi- 
cult to  calculate,  but  are  known  to  exist  since  the  attenu- 
ation does  not  vanish  when  the  fluid  is  removed.  Since  the 
viscous  attenuation  depends  on  the  permeability,  the  frame 
1 o 8 s can  be  expected  to  dominate  in  low  permeability  clays 
and  rocks.  If  the  frame  loss  is  linear  with  frequency, 
then,  to  first  approximation,  so  will  be  the  bulk  attenua- 
tion. Second,  the  theory  is  based  on  a uniform  pore  size. 
If  the  sediment  is  poorly  sorted,  then  the  transition  zone 
between  the  f2  and  f1/2  dependence  may  extend  over  a wide 
band  of  frequencies.  Third,  only  measurements  in  the  same 
sample  can  be  compared  since  permeability  changes  between 
slightly  different  particle  sizes  will  obscure  the  pre- 
dicted effects.  A study  in  a single  sand  over  an  extended 
frequency  range  is  lacking.  However,  as  noted  earlier,  at- 
tenuation in  well  sorted  sand  samples  has  been  observed  to 
vary  as  f1/2.  An  observation  of  the  small  amount  of  dis- 
persion implied  by  the  theory  will  require  accurate  mea- 
surements over  a frequency  band  that  includes  the  transi- 
tion zone.  Most  measurements  have  been  made  at  either  very 
low  or  very  high  frequencies. 
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Biot  (1962)  recognized  that  losses  due  to  the 
frame  must  be  included  in  a general  theory  of  acoustic 
propagation  in  porous  media.  He  suggested  that  the  elastic 
moduli  be  replaced  by  complex  operators  and  outlined  the 
form  of  possible  distributions  of  retardation  times.  Stoll 
(197^  , 1977)  and  Stoll  and  Bryan  (1970)  have  made  extensive 
quantitative  studies  based  on  the  simplifying  assumption 
that  the  operators  can  be  approx*' . 'ted  by  complex  constants 
independent  of  frequency.  In  effect,  this  adds  a constant 
frame  log  decrement  to  the  viscous  loss. 

Several  authors  have  studied  special  aspects  of 
the  general  viscous  interaction  model.  Geertsma  and  Smit 
(1961),  Deresievicz  and  Rice  (1962,  196U),  and  Horton  (1978) 
have  considered  the  effects  of  boundaries.  Brutsaert 

V 

(1963)  has  studied  partially  saturated  systems.  Smith  and 
Greenhorn  (1972,  Smith  et  al . , 197**)  have  considered  ther- 
mal conductivity  as  well  as  viscous  dissipation  in  a rigid 
frame.  Mengi  and  McNiven  (1977,  1978)  have  studied  the 
response  of  porous  media  to  a transient  pulse  and  the  in- 
fluence of  anisotrophy  on  acoustic  wave  shapes.  Chase 
(1979)  has  applied  the  theory  to  saturated  foam. 

The  next  section  presents  a simplified  derivation 
of  the  equations  that  govern  shear  wave  propagation  in 
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in  porous  media-  The  end  result  is  equivalent  to  Stoll's 
(1977)  formulation  of  the  theory. 

Viscous  Loss  Model 


The  displacement  of  a damped  harmonic  oscillation 
is  given  by 

u = ei(u>t-*x)  , III-l 


where  Real  ( 5, ) = £ = ^ “ y 

Img  (T)  = Z ' = a (attenuation  in  nepers/unit 

distance)  . 


The  elastic  wave  equation  for  shear  waves,  that  is,  dis 
placements  perpendicular  tq  the  direction  of  motion,  is 


y 


8 2u 
Z 

3t2 


III-2 


Substitution  of  Eq.  III-l  into  Eq.  III-2  yields 


y*2  = u2p  III-3 

which,  in  the  absence  of  dissipation,  gives  the  familiar 
expression  for  the  shear  wave  speed  in  terms  of  the  shear 
modulus  and  the  density. 


— jjg 
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III-U 


v2  = p/p  . 


However,  if  Eq.  III-3  is  applied  for  i ' ^ 0,  then  either  p 
or  p,  or  both,  must  be  considered  as  complex  operators. 

If  p is  arbitrarily  assigned  an  imaginary  compo- 
nent, then  Eq.  III-3  implies 


III-5 


III-6 


III-7 


where  r = aV/w. 

When  the  damping  is  small,  r2>>l  and  Eqs . III-5  and  III-7 


reduce  to 


V2  = p/p 
6 = aV/f  . 


III-8 


Equation  III-8,  which  is  independent  of  the  frequency,  cor- 
responds to  the  elastic  expression  for  the  speed.  Equa- 
tion III-9  yields  a constant  log  decrement  if  a varies  as 
the  first  power  of  the  frequency.  This  is  the  viscoelastic 
model  favored  by  Hamilton  (1976b).  The  equation  given  in 


Fig.  II-l  are  valid  when  Eq.  III-7  can  be  approximated  by 
Eq.  III-9 . 


The  above  expressions  pertain  to  the  properties 
of  the  bulk  system.  However,  since  water  has  vanishing  ri- 
gidity, Gassman  (1951)  concluded  that 

li(bulk)  = y(frame)  . III-10 

Rigorously,  the  properties  of  the  frame  are  de- 
fined only  in  terms  of  the  total  sediment.  However,  it  is 
often  assumed  in  laboratory  studies  that  values  obtained 
from  the  dry  sediment  can  be  used  to  approximate  the  frame 
parameters  when  the  sediment  is  saturated  (e.g.,  Gardner 
et  al . , 1964;  Gregory,  1976;  Domenico,  1977)*  It  is  then 
possible  to  predict  the  effects  of  fluid  saturation. 

From  small  damping,  application  of  Eq.  III-10  to 
Eqs.  III-8  and  III-9  implies 


where  V(frame)  is  the  shear  speed  in  the  dry  material. 


Equation  III-ll  states  that  the  bulk  shear  speed 
is  less  than  that  of  the  frame  in  proportion  to  the  square 
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root  of  the  ratio  of  the  two  densities.  Since  dissipation 
has  been  introduced  only  through  the  imaginary  part  of  the 
shear  modulus,  which,  for  the  moment,  is  assumed  equal  for 
both  the  frame  and  the  hulk  sediment,  the  attenuation  in 
each  case  is  the  same.  This  leads  to  Eq.  III-12  which 
states  that  the  bulk  log  decrement  is  less  than  that  of  the 
frame . 

Data  in  the  literature,  in  contradiction  to 
Eq.  1 1 1-12 , indicate  that  attenuation  increases  with  satu- 
ration. This  implies  an  additional  loss  mechanism  due  to  the 
presence  of  the  fluid.  Conceivably,  the  fluid  could  influence 
the  frictional  losses  at  the  grain  boundaries.  That  is, 
the  frame  loss  in  the  dry  sediment  could  be  different  from 
the  frame  loss  in  the  saturated  sediment.  Alternatively, 
the  additional  attenuation  could  arise  from  viscous  drag 
caused  by  the  relative  movement  of  the  frame  and  the  pore 
fluid . 

Equation  III-ll  was  obtained  with  the  assumption 
that  the  entire  mass  of  the  fluid  was  perfectly  coupled  to 
the  frame  (Gassman's  "closed  system").  However,  if  the 
coupling  is  due  to  viscous  interaction,  then  the  possibili- 
ty exists  that  only  a portion  of  the  fluid  will  move  with 
the  frame.  The  p(bulk)  in  Eq.  III-ll  should  then  be 


72 


replaced  by  an  "effective  density"  which  in  general  will 
depend  upon  the  frequency.  If  the  coupling  is  less  than 
perfect , then  a smaller  change  between  the  dry  frame  and 
the  bulk  velocities  is  indicated  than  that  predicted  by 
Eq.  III-ll.  This  has  been  observed  by  Hardin  and  Bichart 
(1963)  in  Ottawa  and  crushed  quartz  sand  and  by  Domenico 
(1977)  in  pressurized  sand  and  glass  bead  specimens. 

Ament  (1953)  obtained  an  equation  for  the  effec- 
tive density  based  on  a simple  consideration  of  viscous, 
incompressible  flow.  He  found 

p(effective)  - pUtatic)  * j^j - HI-13 

where  K is  the  absolute  permeability,  ri  the  fluid  vis- 
cosity, B the  porosity,  and  A = pg  - pf  is  the  difference 
between  the  grain  and  fluid  densities. 

The  imaginary  part  of  Eq.  III-13  will  introduce 
an  additional  attenuation  due  to  the  relative  movement  of 
fluid  and  frame  that  at  first  increases  as  the  square  of 
the  frequency  and  then  becomes  almost  constant.  In  order 
for  the  permeability  measured  at  steady  flow  rates  to  be 
appropriate.  Ament  states  that  the  frequency  should  be  less 
than  approximately 
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u>  = 2n/pfa2  , 


where  a is  the  pore  "size." 

In  a detailed  study  of  relative  frame-fluid  move- 
ment caused  by  acoustic  stress,  Biot  (1956b)  obtained  an 
expression  similar  to  Eq.  III-lU.  The  assumption  of  Pois- 
euille  flow  is  invalid  above  this  frequency.  He  derived  a 
complex  correction  factor  for  the  viscosity  that  allows 
consideration  of  higher  frequencies.  His  result,  in  terms 
of  zero  order  Kelvin  functions,  was 

KT(  k ) 


K(„]  - 1 *TU) 

t 1 + 2T  ( k ) i / 1 


III-15 


where 


x _ ber1 (k)  + ibei ' ( < ) 
v ber ( k ) + ibei ( k ) 


III-16 


and 


k = a(u>pf/n) 


l/2 


III-17 


The  pore  size  a is  strictly  the  radius  of  a cir- 
cular tube.  However,  Biot  found  that  the  analytic  solution 
for  flow  between  two  infinite  parallel  plates  gave  essen- 
tially the  same  correction  if  a was  adjusted  by  a scale 
factor.  Thus  a may  be  regarded  as  a pore  size  parameter 
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that  will  vary  with  pore  shape  and  sinuosity.  However,  a 
large  variation  in  pore  size  is  not  allowed  in  the  theoret- 
ical formulation. 

Replacing  the  viscosity  in  Eq.  13  by  nF(tc)  as  de- 
fined in  Eq.  III-2U  results  in  a viscous  loss  term  that  in- 
creases as  the  square  of  the  frequency  at  lower  frequencies 
and  as  the  square  root  of  the  frequency  at  high  frequen- 


cies. 


gives 


With  both  jj  and  p regarded  as  complex,  Eq.  III-3 


*2  = 2- 


2 yz  + y 


1 + /I  + 


( p ' y + py  ' V 
\ py+p  ' y ' / 


III-18 


<2  - zUL 


y ^ + y ' 2 


/ p 1 y + py 1 \2] 
\ py+P ' V ' / I 


. III-19 


If  y is  taken  as  an  empirical  complex  constant 
determined  by  the  frame,  and  p is  Ament's  effective  den- 
sity, Eq.  III-13,  modified  by  Biot's  frequency  correction 
to  the  viscosity,  Eq.  III-15,  then,  with  the  appropriate 
choice  of  parameters,  Eqs.  III-18  and  III-19  are  equivalent 
to  Stoll's  (1978)  application  of  Biot's  theory  for  shear 
wave  propagation  in  porous  media.  This  is  not  unexpected 
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since  Geersma  and  Smit  (1961)  showed  that  Ament's  equation 


| 


I 

\ 


f. 


can  be  derived  from  Biot's  equations  of  motion.  The  only 
difference  between  Eqs.  III-18  and  III-19  and  similar  form- 
ulas used  by  Stoll  is  the  number  of  parameters  which  must 
be  specified.  Stoll's  formulation  includes  not  only  the 
pore  size  parameter  discussed  above,  but  also  an  additional 
"structure  constant." 

The  structure  constant  in  effect  adds  an  addi- 
tional term  to  the  effective  density  of  Eq.  III-13  to  ac- 
count for  the  possibility  of  mass  coupling  between  the 
frame  and  fluid  not  due  to  viscous  interaction,  such  as 
would  exist  in  isolated  cavities  oriented  perpendicular  to 
the  pressure  gradient  (Zwikker  and  Kosten,  19U9).  For  the 
highly  permeable  sediments  considered  in  the  present  study, 

V 

stagnant  pockets  of  fluid  can  be  neglected.  This  amounts 
to  setting  Stoll's  structure  constant  to  unity  or  Biot's 
(1956a)  mass  coupling  parameter,  Pi2»  tc  zero.  Slight 
variations  of  these  parameters  cause  only  minor  changes  in 
the  shape  of  the  resulting  attenuation  versus  frequency 
curves . 

Figure  III-l  illustrates  the  relative  theoretical 
effects  of  the  various  alterations  to  the  elastic  solution 
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of  the  shear  wave  equation  of  motion.  The  values  are 
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appropriate  for  glass  beads  .18  mm  in  diameter.  With 
p = p ( dry ) , a constant  log  decrement  for  the  frame,  given 
by  y',  yields  a straight  line  with  unity  slope.  With 
y'  = 0,  the  viscous  loss  for  constant  flow  resistance  ob- 
tained from  Ament's  effective  density  gives  a constant  at- 
tenuation at  high  frequencies.  Introduction  of  Biot's  fre- 
quency correction  to  the  viscosity  causes  the  viscous  loss 
at  high  frequencies  to  vary  as  the  square  root  of  the  fre- 
quency. The  total  predicted  attenuation  when  all  parame- 
ters in  Eq.  III-19  are  finite  is  the  sum  of  the  individual 
frame  and  viscous  losses. 

Figure  III-2  shows  the  variation  of  the  predicted 
viscous  loss  with  changes  of  grain  size  and  permeability. 
The  input  parameters  correspond  to  expected  values  for  four 
sizes  of  glass  beads  near  the  surface  of  the  sample.  The 
transition  region  between  the  asymptotic  f2  and  fi/2  de- 
pendence of  attenuation  on  frequency  is  seen  to  depend  upon 
the  permeability,  which  in  turn  is  controlled  by  the  grain 
size.  The  transition  occurs  at  lower  frequencies  for  larg- 
er grain  sizes,  in  accord  with  Eq.  Ill-lit.  At  low  frequen- 
cies the  viscous  attenuation  is  expected  to  increase  with 
grain  size.  At  higher  frequencies,  a simple  statement  of 
the  relation  between  grain  size  and  attenuation  is  no 
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FIGURE  I I 1-2 

VARIATION  OF  PREDICTED  VISCOUS  LOSS 
WITH  GRAIN  SIZE  AND  PERMEABILITY 
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longer  possible  since  the  frequency  dependent  breakdown  of 
Poiseuille  flow  causes  the  curves  to  overlap. 

Figure  III-3  illustrates  the  theoretical  depen- 
dence of  shear  wave  speed  on  frequency.  The  input  parame- 
ters correspond  to  the  four  sizes  of  glass  beads  considered 
in  Fig.  III-2.  The  absolute  value  of  the  speed  depends  on 
the  value  assigned  to  the  real  part  of  the  shear  modulus. 
This  value  was  purposely  varied  with  grain  size  in  order  to 
simplify  the  display.  Therefore  the  variation  with  grain 
size  of  the  absolute  magnitude  of  the  phase  velocity  shown 
in  the  figure  is  not  a prediction  of  the  theory.  The  rela- 
tive change  in  speed  with  frequency  is  the  quantity  of  in- 
terest. For  a given  grain  size,  the  speed  at  low  frequen- 
cies corresponds  to  perfect  coupling  between  the  fluid  and 
the  frame.  That  is,  the  density  term  in  the  expression 
V2  = y/p  is  the  total  density  of  the  sediment.  At  infinite 
frequency,  the  fluid  no  longer  moves  with  the  frame.  The 
limiting  value  of  the  speed  is  then  given  by  the  density  of 
the  frame  alone.  The  overall  change  in  speed  between  high 
and  low  frequencies  depends  upon  the  porosity,  that  is,  the 
ratio  of  the  volume  of  the  fluid  to  the  total  volume.  Typ- 
ical values  of  porosity  represent  a predicted  dispersion  of 


FIGURE  III-3 

PREDICTED  SHEAR  RAVE  DISPERSION  FOR  VARIOUS  GRAIN  SIZES 
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a few  percent  from  zero  to  infinite  frequency.  Fig- 
ure III-3  shows,  however,  that  most  of  the  change  occurs  in 
a relatively  narrow  frequency  hand,  and  that  this  region 
depends  upon  the  grain  size  and  permeability. 
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CHAPTER  IV 

SEDIMENT  PHYSICAL  PROPERTIES 

The  review  presented  in  Chapter  II  of  the  shear 
wave  data  available  in  the  literature  reveals  a scarcity  of 
measurements,  particularly  of  attenuation  in  unconsolidated 
porous  media.  There  are  sufficient  data  to  suggest  trends 
in  the  variations  of  shear  wave  parameters  with  changes  in 
sediment  physical  properties,  but  the  empirical  relations 
so  established  are  inadequately  determined.  Much  of  the 
current  knowledge  of  shear  wave  propagation  in  marine  sedi- 
ments is  based  on  measurements  performed  in  partially  sat- 
urated soils  or  in  rocks.  While  these  are  expected  to  par- 
allel approximately  the  acoustic  behavior  of  completely 
saturated  unconsolidated  sediments,  the  possibility  exists 
that  the  mechanisms  responsible  for  attenuation  may  differ 
significantly.  It  is  possible  to  predict  theoretically  the 
effects  of  certain  loss  mechanisms.  However,  the  available 
data  are  inadequate  to  test  the  models  proposed.  Supple- 
mentary data  are  supplied  in  the  present  study. 

In  order  to  monitor  the  relationships  between  an 
acoustic  parameter,  such  as  the  wave  speed,  and  a physical 
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The  Panama  City,  or  PC,  sand  is  a mature,  medium 
white  quartz  beach  sand  from  Panama  City,  Florida.  The  Ot- 
tawa sand  is  a supermature,  medium  to  coarse  quartz  sand 
sold  as  a standard  for  field  tests  of  soil  density.  Tex- 
tural maturity  reflects  the  sorting  and  smoothness  of  the 
grains.  The  Ottawa  grains  are  fairly  well  rounded,  the  PC 
grains  are  more  angular.  The  Red  sand  is  a submature, 
coarse,  riverbed  sand  of  mixed  quartz  and  fledspar  grains. 
Values  of  the  density  and  bulk  modulus  of  pure  quartz  are 
2.65  g/cm3  and  3.8  * 1010  N/m2  (Gregory,  1975)* 

The  viscous  interaction  model  outlined  in  the 
previous  section  requires  specification  of  several  physical 
parameters.  These  are  the  frame  and  fluid  densities,  the 
porosity,  the  permeability,  and  the  fluid  viscosity.  Water 

V 

was  used  as  the  pore  fluid.  It  has  a density  of  1 g/cm3 

p 

and  a viscosity  of  approximately  .001  nsec/m  depending  on 
the  temperature  ( CRC , 1963).  The  densities,  porosities, 
and  permeabilities  of  the  sediments  as  measured  in  the  lab- 
oratory are  given  in  Table  I. 

The  density  and  porosity  of  a regular  packing  of 
identical  spheres  should  not  change  with  grain  diameter. 

The  expected  value  of  porosity  for  a random  packing  of  like 
spheres  is  .38  (Carman,  1956).  The  fluctuation  of  the 
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TABLE  I 


SEDIMENT  PHYSICAL  PROPERTIES 


Sediment 

Grain 
Size 
(mm ) 

Porosity 

Dry 

Density 
( g/cm3 ) 

Saturated 

Density 

Permeability 
(m2  ) 

L 

.0625 

.389  ± .004 

1.42  ± .008 

1.81 

2.58  x 10-12 

MH 

.177 

.365  ± .008 

1.55  ± .02 

1.92 

1.70  x 10_n 

XHX 

.250 

.369  ± .007 

1.55  ± .01 

1.92 

—3 

ON 

X 

H 

O 

t— • 

XPX 

• 707 

.373  ± .006 

1.53  ± .009 

1.90 

4.63  X 10_1° 

RED 

— 

.383 

1.64 

2.02 

5.71  X 10"u 

PC 

• 35^ 

.386  ± .014 

1.60  ± .03 

1.99 

3.79  x 10"11 

OTT 

• 50 

3.56  ± .005 

1.69  ± .01 

2.05 

1.39  x 10“19 

87 


measured  values  about  the  ideal  porosity  reflect  variations 
in  grain  size  and  shape.  However,  to  first  order,  the  po- 
rosities are  the  same.  The  densities  and  porosities  were 
determined  by  adding  a known  volume  of  water  to  a weighted 
volume  of  sediment  packed  by  vibration.  The  standard  de- 
viations of  four  separate  measurements  are  given  in  the  ta- 
ble. 


A constant  head  permeameter  was  constructed  ac- 
cording to  ASTM  standards  (ASTM  D2l*3^-68)  to  obtain  values 
of  the  permeability.  Figure  IV-2  is  a scale  drawing  of  the 
apparatus.  A differential  pressure  gauge  was  used  to  mea- 
sure the  pressure  drop,  h,  due  to  the  flow  resistance  of 
the  sediment.  The  coefficient  of  permeability,  k,  is  re- 
lated to  the  quantity  of  water,  Q,  that  flows  in  time,  t, 

■by 


k = QL/thA 


IV-1 


where  L is  the  sample  length  and  A the  cross-sectional  ar- 
ea. Figure  IV-3  is  a sample  plot  of  flow  rate,  Q/At,  ver- 
sus the  pressure  gradient,  h/L.  The  slope  of  the  line  is 
the  permeability,  k,  in  cm/sec  at  the  temperature  T of  the 
experiment.  The  ratio  of  the  viscosity  of  water  at  T to 
that  at  20°C  was  used  to  reduce  the  measurements  to  a stan- 
dard temperature-viscosity  reference. 


FIGURE  IV- 3 

DETERMINATION  OF  SEDIMENT  PERMEABILITY 
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The  influence  of  fluid  viscosity  and  density  can 
he  removed  from  the  permeability  constant  k by  defining  an 


absolute  permeability 


K = kn/pg  • 


IV-2 


where  tj  is  the  fluid  viscosity,  p the  density,  and  g the 

acceleration  of  gravity.  For  water  at  20°C,  K(m2)  = 

_ 3 

1.02  * 10  k(m/sec).  Measured  values  of  the  absolute  per- 
meability in  m2  are  given  in  Table  I. 

The  data  in  Table  I reveal  that  permeability  in- 
creases with  increasing  grain  size.  This  is  a general  ob- 
servation for  flow  through  unconsolidated  materials  (Car- 
man, 1956).  The  permeability  can  be  expressed  in  terms  of 
the  square  of  the  mean  grain  diameter  by  the  Kozeny-Carman 
equation  for  uniform  spheres 


„ . si. 

(1_s,2 


IV-3 


where  8 is  the  sediment  porosity.  The  factor  of  36  is  ob- 
tained from  the  specific  surface  area  of  spherical  grains 
and  is  sometimes  assigned  different  values  for  other  parti- 
cle shapes.  The  empirical  parameter  kQ  varies  from  4.5 
5.1  for  gas  flow  through  packings  of  spheres  (Carman, 

1956).  The  value  obtained  from  the  present  data  is  4.0  for 
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the  heads  and  5.6  for  well  sorted  sands  (Fig.  IV-k ) . Addi- 
tional sand  data  were  obtained  from  Nolle  et  al.  (1963). 

Both  the  specific  surface  area,  represented  by 
the  value  36,  and  the  parameter  kQ  are  necessary  to  adapt 
equations  established  for  flow  in  parallel  tubes  of  con- 
stant cross-section  to  the  irregular  pore  shapes  found  in  a 
loose  aggregate.  The  flow  resistance  depends  on  the  sur- 
face area  exposed  to  the  fluid  and  the  rate  at  which  the 
fluid  moves  past  the  surface.  The  flow  rate  measured  in  a 
permeability  test  is  actually  less  than  the  fluid  speed  in 
the  medium  due  to  the  sinuosity  of  the  pores.  That  is,  the 
path  followed  by  the  fluid  is  longer  than  the  length  of  the 
sample.  A constant  value  of  kQ  for  the  various  sediments 
under  consideration  implies  that  the  effects  of  sinuosity 
can  be  characterized  by  a single  number. 

Biot  (1956)  obtained  a high  frequency  correction 
to  the  viscous  resistance  in  a porous  media  induced  by 
acoustic  stress  by  consideration  of  the  velocity  pattern 
within  a uniform  circular  tube.  His  formula  included  the 
tube  diameter  (Eq.  III-17)-  It  was  suggested  that  the  val- 
ue of  the  diameter  could  be  empirically  adjusted  to  account 
for  the  effects  of  pore  shape  and  sinuosity.  The  Kozeny- 
. Carman  equation  allows  Biot’s  pore  size  parameter  to  be 
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expressed  in  terms  of  the  mean  grain  diameter  determined  by 
permeability  measurements  at  low  flow  rates.  Hovem  (1979) 
concluded  that 


a 


d B 
3 (1-B) 


IV-U 


where  a is  the  pore  size  parameter  in  Eq.  III-17  and  d is 
the  mean  particle  size  in  Eq.  IV-3.  This  compares  favor- 
ably with  Stoll's  (1976)  suggested  value  of  1/6  the  parti- 
cle size. 

The  usefulness  of  Eq.  IV-U  depends  upon  the  va- 
lidity of  the  Kozeny-Carman  equation.  For  certain  materi- 
als, including  many  rocks  (Carman,  1956),  kQ  is  not  a con- 
stant. That  is,  the  effects  of  pore  shape  and  sinuosity 
are  not  easily  related  to  a grain  or  pore  size.  For  un- 

V 

consolidated  sands,  however,  Eq.  IV-U  provides  a value  of 
the  pore  size  parameter  independent  of  acoustic  measure- 
ments . 


L 


i __ 


i 


C H A P T E R V 


TRANSDUCER  CHARACTERISTICS 
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A primary  reason  for  the  scarcity  of  shear  wave 
measurements  in  marine  sediments  has  been  the  lack  of  a 
source-receiver  combination  capable  of  detecting  such  waves 
in  low  rigidity  media.  The  use  of  bender  elements  has 
overcome  this  difficulty  (Shirley  and  Anderson,  1975a; 
Shirley  and  Hampton,  1978).  A bender  element  consists  of 
two  ceramic  plates  cemented  back  to  back  but  oppositely  po- 
larized. An  applied  voltage  causes  one  plate  to  contract 
while  the  other  expands.  This  causes  the  assembly  to  bend 
in  a manner  similar  to  the  action  of  a bimetal  strip  in  a 
thermostat  (Fig.  V-l).  Displacement  of  the  sediment  per- 
pendicular to  the  length  of  the  element  causes  a shear  wave 
to  propagate  parallel  to  the  axis  of  the  transducer.  Bend- 
er elements  are  commercially  available. 

Previous  laboratory  work  (Shirley  and  Anderson, 
1975b;  Shirley,  1977)  has  been  conducted  with  transducer 
arrays  of  from  one  to  six  bender  elements  separated  by  com- 
pliant spacers.  Only  the  ends  of  the  elements  were  direct- 


ly coupled  to  the  sediment 


The  present  research  has 

9k 


DIRECTION  OF 

FLECTRICAL 

POLARIZATION 


FIGURE  V-1 

OPERATION  OF  A CERAMIC  BENDER  TRANSDUCER 

(after  Shirley  1978) 
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utilized  an  improved  laboratory  design  (Shirley,  1978)  in 
which  the  entire  length  of  the  element  drives  the  sediment. 
This  significantly  increases  the  signal  amplitude  for  a 
given  flexion  of  the  transducer.  The  bender  elements  are 
suspended  beneath  a thin  stainless  steel  plate  by  a semi- 
flexible  urethane  potting  compound.  The  casting  material 
provides  both  support  and  electrical  insulation,  but  is 
pliable  enough  to  allow  the  elements  to  vibrate  freely. 
Figure  V-2  shows  the  two  transducer  mounts  used  in  the 
present  study. 

The  free  air  displacement  at  the  end  of  a bender 
element  is  given  by  (Shirley  and  Anderson,  1975a) 

D = 1.5d31L2V/T2  , V-l 

where  L and  T are  the  length  and  thickness  of  the  element, 

V is  the  applied  voltage  and  d3i  is  the  piezoelectric  con- 
stant of  the  ceramic  material.  The  individual  ceramic 
plates  are  operated  in  the  length  extensional  mode,  that 
is,  the  applied  voltage  is  across  the  thickness  of  the 
plate  while  the  resultant  mechanical  strain  is  parallel  to 
the  length  of  the  plate.  The  subscripts  3 and  1 simply  in- 
dicate that  the  strain  is  perpendicular  to  the  direction  of 
polarization. 


■ 
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The  bender  elements  used  in  the  present  study 


have  a free  air  displacement  of  approximately  10~4  cm.  The 
actual  displacement  should  be  less  when  the  element  is 
loaded  by  the  potting  compound  and  the  sediment.  The  dis- 
placement is  sufficient  to  generate  shear  waves  in  low  ri- 
gidity sediments  but  is  small  enough  to  give  a linear  re- 
sponse. A variation  of  acoustic  parameters  with  displace- 
ment amplitude  was  not  observed.  Hardin  and  Richart  (1963) 
noticed  amplitude  effect  for  torsional  vibrations  between 
10~4  and  10~2  radians  of  angular  displacement.  Figure  V-3 
Illustrates  the  quality  of  a received  shear  wave  pulse. 

The  bender  elements  are  low  Q devices  with  a free 
air  resonance  that  is  inversely  proportional  to  the  square 
of  the  length  (Shirley  and  Anderson,  1975a).  However,  the 
shear  wave  speed  in  the  sediment  controls  the  frequency  re- 
sponse vhen  the  element  is  loaded. 

A knowledge  of  the  transducer  frequency  response 
would  be  valuable  in  the  determination  of  attenuation  as  a 
function  of  frequency.  An  accurate  characterization  re- 


quires the  transducer  to  be  mounted  in  the  sediment  since 
the  response  depends  upon  the  load.  This  in  turn  requires 
the  frequency  response  of  the  medium  to  be  known.  However, 
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the  response  of  the  sediment  is  the  goal  of  the  present 
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re- 
search. There  has  not  yet  been  sufficient  work  done  to  es- 
tablish a reference  medium  whose  shear  wave  properties  are 
known  accurately  enough  to  allow  routine  calibration  of  the 
elements.  It  is  possible,  however,  to  solve  for  the  trans- 
ducer frequency  response  once  the  sediment  attenuation  has 
been  determined. 

Figure  V-U  shows  the  measured  amplitude  response 
of  a 2.5^  cm  bender  element  in  three  different  sediment 
types.  The  attenuation  in  each  sediment  was  first  deter- 
mined from  the  amplitude  decay  with  transducer  separation 
at  fixed  frequencies.  The  measured  amplitudes  at  a single 
separation  were  then  corrected  to  give  the  response  of  the 
two-element  projector-receiver  combination.  Figure  V-5 

V 

shows  the  variation  with  transducer  separation  of  the  ob- 
served output  voltages  in  the  PC  sand. 

Figure  V-U  clearly  shows  the  dependence  of  the 
transducer  response  on  the  wavelength  excited  in  the  sedi- 
ment. The  measured  shear  wave  speed  decreases  from  sand  to 
clay.  The  figure  also  indicates  the  useful  bandwidth  of 
the  transducers.  The  small  bender  elements  generate  mea- 
surable signals  over  approximately  a decade  of  frequency. 
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There  are,  however,  several  practical  considerations  which 
limit  the  useful  range  to  less  than  a decade. 

A reliable  determination  of  the  speed  and  the  at- 
tenuation requires  measurements  at  several  transducer  sepa- 
rations. However,  the  increase  of  sediment  attenuation 
with  frequency  shifts  the  maximum  received  amplitude  lower 
in  frequency  as  the  separation  of  the  transducers  is  in- 
creased. Thus  the  lower  frequencies  of  the  transducer  out- 
put are  selectively  favored  even  though  more  energy  may  be 
transmitted  at  a higher  frequency.  The  effect  is  easily 
seen  in  Figure  V-5.  At  1 cm,  the  maximum  output  voltage 
occurs  at  Just  over  4 kHz,  which  is  close  to  the  resonance 
of  the  transducer  in  sand.  However,  at  8 cm  separation, 
the  maximum  amplitude  occurs  near  3 kHz.  In  sediments  with 

V 

high  attenuation,  the  upper  frequencies  are  lost  in  the 
noise  level  before  a sufficient  number  of  readings  at 
increased  separations  can  be  made  with  which  to  calculate 
confidently  the  speed  and  attenuation. 

The  finite  size  of  the  test  tank  limits  the  use- 


fulness  of  the  lower  frequencies  of  the  transducer  band- 
width. The  direct  pulse  between  the  transducers  must  be 
relatively  short  in  order  to  distinguish  the  direct  arrival 


shows  a signal  returned  from  the  bottom  of  a 16  x 20 
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30  cm  aluminum  tank  that  arrives  Just  after  the  direct 
event.  At  lower  frequencies  the  two  will  overlap.  The 
driving  signal  must  be  gated  to  produce  a short  pulse.  In 
the  photograph,  only  a single  cycle  of  the  specified  fre- 
quency is  input  to  the  transducer. 

Measurements  on  a pulsed  signal  do  not  refer 
strictly  to  a single  frequency;  rather  they  apply  to  a band 
of  frequencies  controlled  by  the  width  of  the  pulse  and  the 
characteristics  of  the  filtering.  If  the  transducer  re- 
sponse sharply  increases  over  the  bandwidth  of  the  filtered 
signal,  then  the  received  signal  amplitude  will  be  domi- 
nated by  the  higher  frequencies  in  the  band  rather  than  by 

V 

the  center  frequency  of  the  input  signal.  This  effect  can 
be  observed  on  the  oscilloscope  trace.  In  Figure  V-3  the 
received  signal  has  essentially  the  same  period  as  the  sin- 
gle cycle  driving  pulse.  However,  at  low  frequencies  the 
period  of  the  received  signal  is  noticeably  less  than  that 
of  a continuous  signal  at  the  assigned  frequency.  Data 
taken  at  the  edges  of  the  transducer  response  are  therefore 
somewhat  ambiguous,  even  though  the  amplitude  is  sufficient 
for  a reading.  Better  discrimination  of  the  signal 
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frequency  can  be  obtained  by  adjustment  of  the  output  fil- 
ter to  a higher  Q.  This,  however,  also  greatly  reduces  the 
received  signal  amplitude. 

An  additional  problem  arises  at  low  frequencies 
due  to  the  unknown  radiation  pattern  of  the  bender  ele- 
ments. The  divergence  of  the  transmitted  energy  due  to 
geometrical  spreading  must  be  removed  from  the  observed 
amplitude  decay  before  a reliable  determination  of  the  sed- 
iment attenuation  can  be  made.  If  the  geometrical  spread- 
ing is  correctly  removed,  then  a plot  of  relative  signal 
amplitude  expressed  in  dB  will  vary  linearly  with  separa- 
tion. The  slope  is  the  attenuation  per  unit  distance  due 
to  the  sediment.  Figure  V-6  shows  the  relative  signal 
amplitudes  obtained  from  Figure  V-5  after  removal  of  as- 
sumed  spherical  spreading  using  the  center  to  center  dis- 
tance between  the  transducers.  This  results  in  a consis- 
tent determination  of  the  attenuation  which  in  turn  allows 
calculation  of  the  transducer  response  shown  in  Figure  V-U. 
However,  an  easily  observed  signal  was  obtained  at  fre- 
quencies less  than  3 kHz.  These  data  are  shown  in  Fig- 
ure V-7  with  the  same  spherical  divergence  correction  ap- 
plied as  in  Figure  V-6.  The  progressive  curvature  of  the 
amplitude  versus  separation  lines  suggests  that  the 
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assumption  of  spherical  spreading  is  no  longer  valid.  The 
curves  become  linear  at  greater  and  greater  separations  as 
the  wavelengths  increase.  The  wavelengths  associated  with 
the  frequencies  in  Figure  V-7  range  from  5 to  20  cm.  That 
is,  for  low  freqeuncies  the  separations  allowed  in  the 
small  tank  are  much  less  than  a wavelength  from  a transdu- 
cer which  itself  has  a definite  length.  The  implication  is 
that  the  nearfield  of  the  transducer  has  been  encountered. 
Consistent  measurements  at  lower  frequencies  in  the  silt 
and  the  clay  are  presumably  possible  because  of  the  lower 
speeds  and  hence  shorter  wavelengths. 

The  above  considerations  limit  reliable  shear 
wave  measurements  with  the  small  bender  elements  to  slight- 
ly less  than  a decade  of  frequency.  Even  though  such  a 
range  compares  favorably  with  the  bandwidth  of  other  in- 
vestigators, e.g..  Kudo  and  Shima  (1970)  in  Fig.  II-3,  a 
broader  range  was  desired  to  test  the  validity  of  the  vis- 
cous model  of  attenuation.  A study  of  the  preliminary  in- 
vestigations reviewed  in  this  section  indicates  two  means 
of  increasing  the  available  bandwidth.  First,  a larger 
tank  would  permit  both  a longer  driving  pulse  and  greater 
separations.  Second,  since  the  resonance  frequency  of  the 
transducer  depends  on  the  length  of  the  element  , an  array 
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of  different  sized  elements  could  extend  the  frequency 
range.  Both  of  these  improvements  have  been  implemented. 

Figure  V-2  shows  a transducer  that  consists  of 
both  a 2.54  cm  and  a 5.08  cm  bender  element.  These  trans- 
ducers were  used  in  an  aluminum  tank  46  x 56  x 71  cm.  The 
frequency  of  greatest  sensitivity  is  slightly  lower  in  the 
larger  element  than  in  the  shorter  one.  However,  the  Q of 
the  larger  element  is  such  that  its  usable  bandwidth  en- 
compasses that  of  the  smaller  element.  Measurements  with 
the  larger  elements  can  be  obtained  over  slightly  more  than 
a decade  of  frequency.  The  smaller  elements  serve  as  a 
consistency  check. 
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CHAPTER  VI 

LABORATORY  EQUIPMENT  AND  PROCEDURE 

Shear  and  compressional  wave  speed  and  attenua- 
tion were  measured  for  each  of  the  seven  sediments  de- 
scribed in  Chapter  IV,  both  when  dry  and  completely  water 
saturated.  The  transducer  probes  shown  in  Figure  V-2  were 
inserted  in  the  sediments  at  fixed  depth  and  variable  sepa- 
ration. The  projecting  transducers  were  driven  by  a pulsed 
sine  wave  from  one  to  forty  cycles  in  duration.  The  re- 
ceived signal  was  filtered  and  displayed  on  an  oscillo- 
scope screen.  The  change  in  arrival  time  with  separation 
of  a distinctive  peak  of  the  received  signal,  as  measured 

V 

by  a 5-digit  LED  display  function  of  the  oscilloscope,  was 
used  to  compute  the  wave  speed.  Attenuation  was  calcualted 
from  the  decay  of  signal  amplitude  with  separation,  as  mea- 
sured peak  to  trough  on  the  oscilloscope  screen  for  the 
shorter  signals  or  rms  with  a sampling  voltmeter  for  the 
longer  signals.  The  filter  was  adjusted  for  a constant 
Q = 10  to  ensure  that  the  cycle  frequency  of  the  driving 
pulse  was  representative  of  the  received  signal  frequency. 
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A constant  Q,  that  is,  f/Af,  rather  than  a constant  Af  pre- 
serves the  wave  shape  at  different  center  frequencies. 

Measurements  of  the  saturated  sediment  properties 
were  performed  in  a small  aluminum  tank,  16  x 20  x 30  cm, 
using  the  projector-receiver  configuration  shown  in  Fig- 
ure VI-1.  One  of  the  transducer  mounts  was  attached  to  a 
rigid  crossbar  and  remained  in  a set  location.  The  other 
mount  was  fixed  to  a threaded  rod  located  on  top  of  the 
tank.  Rotation  of  this  rod  moved  the  second  set  of  trans- 
ducers relative  to  the  first.  Small  sheets  of  cork  placed 
above  the  transducer  mounts  were  found  necessary  to  prevent 
acoustic  transmission  through  the  support  mechanism.  The 
shear  wave  elements  were  approximately  9 cm  below  the 
sediment  surface. 


Careful  preparation  of  the  sediment  was  necessary 
to  ensure  complete  water  saturation.  The  sands  could  be 
successfully  deaerated  by  slowly  adding  the  sediment  to  the 
water  and  then  vibrating  the  tank  under  vacuum.  The  beads, 
however,  showed  a greater  affinity  for  entraining  air. 

They  were  successfully  deaerated  by  heating  the  sediment  in 


the  tank  to  boiling. 

The  compressional  wave  transducers  mounted  above 
the  shear  wave  elements  allowed  a quantitative  test  of 
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FIGURE  VI-1 

PROJECTOR  RECEIVER  CONFIGURATION 
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complete  saturation.  Only  a small  amount  of  air  trapped  in 
the  sediment  would  significantly  reduce  the  transit  time 
and  amplitude  of  the  compressional  wave  pulse. 

The  compressional  wave  elements  were  thin  piezo- 
electric ceramics  operated  in  the  length  extensional  mode, 
highly  resonant  at  120  kHz.  They  were  driven  by  a two- 
cycle  pulse.  A 100-140  kHz  bandpass  filter  was  applied  to 
the  compressional  wave  output. 

Two  modifications  were  necessary  before  reliable 
shear  wave  results  could  be  obtained  in  the  dry  materials. 
First,  the  transducer  mounts  were  coated  with  conducting 
paint  to  eliminate  electrical  feedover  between  the  source 
and  receiver  circuits  that  obscured  the  beginning  of  the 
shear  wave  pulse.  Second,  a slightly  larger  tank  was  con- 
structed  to  allow  better  discrimination  between  the  direct 
arrival  and  reflection  returns.  This  was  a wood  box  lined 
with  porous  foam  with  inside  dimensions  of  25  x 36  x 24  cm. 
The  same  transducers  and  support  mechanism  were  used  for 
both  the  small  aluminum  tank  and  the  wooden  container.  The 
depth  of  the  transducer  beneath  the  surface  was  the  same  in 
each . 

The  small  tanks  were  used  primarily  to  study  the 
effects  of  grain  size  and  shape  on  the  sediment  acoustic 
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properties.  They  were  convenient  for  this  purpose  since  it 
was  easier  to  saturate  the  smaller  sediment  samples.  Shear 
wave  data  were  collected  at  only  a few  closely  spaced  fre- 
quencies near  the  transducer  resonance.  A larger  tank  was 
used  to  investigate  the  variation  of  shear  wave  parameters 
with  frequency. 

The  large  tank  and  second  set  of  transducers  are 
shown  in  Figure  VI-2.  Desired  use  of  the  tank  as  its  own 
vacuum  chamber  required  a wall  thickness  of  3/1*  in.  (1.9  cm) 
The  tank  volume  was  approximately  150  liters. 

Three  of  the  double  element  transducers  shown  in 
Figure  V-2  were  mounted  in  line  in  the  large  tank.  The 
middle  transducer  was  used  as  the  source,  the  two  outside 
transducers  as  receivers.  The  transducers  could  be  held 
rigidly  at  any  desired  spacing  within  the  tank.  The  depth 
was  set  at  approximately  20  cm  below  the  sediment  surface. 

A refraction  from  the  bottom  of  the  unlined  tank 
restricted  the  use  of  a multicycle  driving  pulse.  A lining 
of  3 cm  thick  porous  foam  with  randomly  spaced  foam  pyra- 
mids 6 cm  in  height  effectively  produced  an  anechoic  cham- 
ber. The  result  is  shown  in  Figure  VI-3.  Noise  interfer- 
ence was  still  present  to  a noticeable  extent  for  the  large 


shear  elements  and  both  elements  showed  some  pulse 
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distortion  at  short  wavelengths.  Most  of  the  measurements 
reported  herein  were  obtained  with  short  driving  pulses. 

The  longer  multicycle  pulses  were  used  primarily  for  qual- 
ity estimation. 

The  variability  of  shear  wave  parameters  with 
only  small  changes  in  sediment  packing  was  of  major  experi- 
mental concern.  Preliminary  investigations  (Shirley  and 
Bell,  1978)  showed  that  self-supporting  structures  were 
possible  in  the  PC  sand  for  porosities  from  roughly  .36 
to  .40.  The  shear  wave  speed  varied  from  130  to  70  m/sec 
with  changes  in  porosity.  The  high  porosity  states  were 
not  stable  but  settled  slowly  with  time  and  rapidly  with 
mechanical  shock.  Since  data  were  taken  at  several  trans- 
ducer separations,  a means  of  returning  the  sediment  to  a 

V 

previous  porosity  was  desired. 

Considerable  research  on  soil  compaction  has  been 
conducted  in  association  with  foundation  design  (Barkan, 
1962;  Wu,  1971)-  Minimum  porosity  packing  of  dry  cohesion- 
less particles  can  be  obtained  by  slowly  pouring  the  mate- 
rial from  a height  into  a container.  The  flow  rate  must  be 
slow  enough  that  each  particle  strikes  the  sediment  sur- 
face independently.  If  such  a packing  is  later  subjected 
to  a small  amplitude  mechanical  shock,  the  porosity  may 
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Increase  slightly.  However,  vibrations  characterized  by 
large  accelerations  will  return  the  dry  sediment  to  minimum 
porosity.  Small  external  stresses  have  little  effect  on 
the  porosity. 

Each  of  the  three  tanks  used  for  shear  wave  mea- 
surements had  vibrators  attached.  The  small  tanks  had  ad- 
justable electromechanical  vibrators  bolted  directly  to 
their  sides.  The  large  tank  was  vibrated  by  an  eccen- 
trically loaded,  single  speed  motor  mounted  beneath  the 
tank  stand.  A noticeable  change  in  both  dry  and  saturated 
sediment  volume  was  observed  upon  initial  vibration.  Ex- 
tended vibration  caused  no  further  reduction  in  volume. 

The  tanks  were  vibrated  during  and  after  a change  in  trans- 
ducer separation.  The  vibration  sufficiently  fluidized  the 

V 

sediment  that  the  thin  transducers  could  easily  be  inserted 
or  moved. 

The  shear  wave  parameters  of  the  sediments  were 
found  to  depend  not  only  on  the  porosity,  but  also  on  small 
internal  stresses  developed  by  the  vibration  process. 

Abrupt  halt  of  the  vibration  can  cause  abnormal  stresses 
within  a loose  aggregate,  particularly  at  the  boundaries  of 
the  container.  Shear  wave  measurements  taken  just  after 
vibration  of  the  tanks  had  longer  delay  times  and  larger 
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amplitudes  than  measurements  obtained  after  the  sediment 
was  allowed  to  relax.  The  decay  in  transit  time  was  ap- 
proximately exponential  and  reached  steady  values  more 
quickly  if  the  vibration  amplitude  was  slowly  decreased 
rather  than  abruptly  stopped.  Compressional  wave  results 
did  not  vary  with  time. 

Measurements  in  the  small  tanks  were  made  after  a 
uniform  settling  time  of  one  and  a half  hours,  after  fif- 
teen minutes  of  low  amplitude  vibration.  The  transducers 
were  active  the  last  fifteen  minutes  of  the  waiting  period. 
Data  in  the  large  tank  were  taken  2h  hours  after  the  trans- 
ducer spacings  were  changed. 

Noticeable  changes  in  the  shear  wave  amplitude 
and  transit  time  could  be  caused  by  very  minor  disturbances 

V 

of  the  sediment,  particularly  for  the  dry  glass  beads  in 
the  small  tank.  Gently  tapping  the  side  of  the  tank  would 
cause  the  transit  time  to  increase.  If  the  disturbance  was 
small,  recovery  of  the  waveform  was  rapid  enough  to  be  vis- 
ible. If  the  disturbance  was  relatively  large,  the  wave 
form  was  distorted  and  a repeat  of  the  vibration  cycle  was 
necessary  to  recover  the  previous  readings.  This  extreme 
sensitivity  presented  the  greatest  obstacle  to  reliable 
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shear  wave  measurements.  The  compressional  wave  signals 
were  unaffected  by  disturbances  of  this  magnitude. 

The  small  compressional  wave  transducers  used  in 
the  saturated  sediments  did  not  couple  sufficiently  well 
with  the  dry  material  for  a compressional  wave  to  be  ob- 
served. To  obtain  estimates  of  the  compressional  wave 
speed  and  attenuation  in  the  dry  sediments,  a small  press 
was  constructed  from  an  aluminum  cylinder  15  cm  in  diameter 
and  30  cm  high.  A projector  and  two  receivers  with  surface 
area  of  approximately  a square  centimeter  were  mounted  on 
short  rods  at  the  bottom  of  the  press  such  that  the  propa- 
gation path  was  across  the  diameter  of  the  cylinder.  The 
two  receivers  were  4.26  and  4.97  cm  from  the  projector  as 
determined  by  calibration  in  three  different  liquids  with 

V 

known  sound  speed.  The  relative  amplitude  difference  be- 
tween the  received  signals  was  also  measured  in  the  fluids. 
Additionally,  a pressure  sensitive  load  cell  was  calibrated 
and  fixed  to  the  bottom  of  the  cylinder.  A signal  of  good 
quality  was  obtained  for  all  of  the  sediments  with  only  a 
small  amount  of  external  frame  pressure.  It  was  hoped  that 
results  for  speed  and  attenuation  could  be  extrapolated  to 
zero  excess  pressure.  However,  the  sensitivity  was  such 


that  only  order  of  magnitude  estimates  could  be  obtained  in 
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that  manner.  Only  values  for  those  sediments  in  which  a 
signal  was  detected  without  compression  are  included  here- 
in. Although  the  received  signals  were  of  poor  quality, 
the  results  were  repeatable  to  within  10  percent  of  one 
another.  A filtered  pulse  was  used  with  a frequency  of  ap- 
proximately 20  kHz. 
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CHAPTER  VII 

DATA  PRESENTATION  AND  ANALYSIS 


Compressional  and  shear  wave  speeds  were  deter- 
mined from  a plot  of  time  delay  versus  relative  transducer 
separation.  Data  at  several  different  spacings  reveal  the 
consistency  of  the  measurement  technique  and  eliminate  the 
need  for  determining  the  exact  beginning  of  the  wave  form 
or  the  absolute  separation  between  the  effective  radiation 
points  of  the  transducers.  Figures  VII-1,  VII-2,  and 
VII-3  present  the  velocity  data  taken  in  the  small  tank. 

The  quoted  values  of  speed  are  obtained  from  a least 
squares  fit  of  the  data.  The  correlations  are  excellent 
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for  the  compressioanl  wave  speeds  but  slightly  less  so  for 
the  shear  wave  speeds  in  the  dry  materials.  Not  shown  are 
points  obtained  in  the  dry  sediments  which  were  widely  di- 
vergent from  the  mean.  In  such  cases  the  vibration  cycle 
was  repeated.  Shear  wave  speed  determinations  were  made  at 
several  frequencies;  only  the  center  frequency  data,  be- 
tween 2 and  4 kHz,  are  shown  in  the  figures.  The  speeds  at 
different  frequencies  differed  by  less  than  5 m/sec  in  the 
dry  sediments. 
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compres sional  wave  attenuation.  The  relative  amplitude  de- 


cay with  transducer  separation  was  first  obtained  in  water. 
Since  water  has  a negligible  attenuation  (-1.3  dB/km  at 
100  kHz;  Urick,  1975  ) this  decay  reflects  the  geometrical 
spreading  pattern  of  the  transducers.  The  correction  fac- 
tor determined  by  this  means  was  applied  to  the  amplitude 
decay  observed  in  the  sediment.  The  result  is  a straight 
line  whose  slope  is  the  attenuation  due  to  the  sediment. 
This  technique  assumes  that  the  transducer  beam  pattern  is 
unaffected  by  the  change  in  impedance  between  water  and 
sediment.  A consistent  linear  decrease  in  the  corrected 
amplitude  with  separation  justifies  this  assumption.  The 
data  shown  in  Figure  VII-U  were  obtained  from  compres s ional 
wave  elements  included  on  a set  of  experimental  shear  wave 


probes.  The  same  technique  was  used  with  the  compressional 
wave  transducers  in  the  small  tank. 

An  analogous  technique  could  be  used  to  calculate 
the  shear  wave  attenuation.  However,  since  the  transducer 
response  varies  with  the  sediment  load,  experimental  deter- 
mination of  the  spreading  loss  separate  from  the  sediment 
attenuation  was  not  feasible.  Spherical  spreading  was  as- 
sumed to  apply  to  the  center  to  center  transducer 
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FIGURE  VI 1-4 

DETERMINATION  OF  COMPRESSIONAL  WAVE  ATTENUATION 


128 

separation.  This  assumption  was  discussed  in  Chapter  7. 
Figure  V-6  shows  the  corrected  attenuation  data  for  the  PC 
sand.  A consistent  interpretation  of  the  data  is  possible. 
However,  Figure  V-6  also  shows  several  values  which  do  not 
lie  on  the  selected  attenuation  curves  even  though  the  time 
delays  associated  with  these  points  were  in  agreement  with 
the  speed  determination.  Such  spurious  points  were  common- 
ly observed  and  hindered  the  interpretation  of  the  shear 
wave  attenuation  data,  particularly  in  the  less  angular 
sediments  which  were  extremely  sensitive  to  minor  distur- 
bances . 

The  first  graph  of  Figure  VII-5  shows  the  rela- 
tive amplitudes  obtained  at  a single  frequency  in  a dry 
sample  of  the  larger  glass  beads.  A spherical  correction 

V 

has  been  applied.  A linear  regression  is  applicable  to  the 
data,  but  the  resulting  correlations  would  be  poor.  Re- 
peated measurements  at  separations  of  2,  3,  and  U cm  show 
more  than  minor  variations  in  the  observed  amplitudes. 
Conceivably,  the  differences  could  arise  from  a change  in 
sediment  attenuation  that  follows  repositioning  of  the 
transducers.  However,  only  points  with  consistent  transit 
times  and  wave  shapes  were  considered.  It  is  unlikely  that 
the  speed  would  remain  constant  while  the  attenuation 


I 


DETE 


130 


varied  "by  a factor  of  2 or  more.  The  assumption  was 
therefore  made  that  the  relatively  large  surface  area  of 
the  tender  elements  gave  rise  to  a coupling  factor  depen- 
dent upon  the  exact  packing  of  the  grains  about  the  trans- 
ducer. This  factor  could  then  be  expected  to  change  when- 
ever the  transducer  was  moved. 

The  received  amplitude  A at  a particular  frequen- 
cy can  be  written  as  the  sum  in  dB  of  the  transducer  fre- 
quency response  T,  the  spreading  loss  SL,  the  attenuation 
times  the  separation  ax,  and  the  hypothetical  coupling 
term  C that  depends  on  the  placement.  That  is, 

A(f,  x)  = T ( f ) + SL(x)  + a ( f ) x + C(x)  . VII-1 

If  the  amplitude  at  a different  frequency  but  the 
same  separation  is  subtracted  from  Eq.  VII-1,  then  an  ex- 
pression is  obtained  for  the  change  in  attenuation  with 
frequency . 

A ( f i ) - A(f2)  • = AT  + Aax  • VII-2 

A plot  of  AA  versus  separation  should  yield  a straight  line 
with  a slope  of  Aa.  The  result  is  shown  in  Figure  VII-5 
for  the  same  data  plotted  on  the  first  graph.  This  sup- 
ports the  assumption  of  a change  in  amplitude  with 
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separation  that  is  independent  of  the  sediment  attenuation 
and  the  frequency. 

The  determination  of  Aa  from  Eq.  VII-2  yields  a 
less  equivocal  interpretation  of  the  amplitude  data  and  al- 
so eliminates  the  necessity  of  assuming  spherical  spread- 
ing. However,  only  a change  in  attenuation,  rather  than  an 
absolute  value,  is  obtained.  Quoted  values  of  attenuation 
are  based  on  the  assumption  that  for  closely  spaced  fre- 
quencies the  attenuation  is  proportional  to  the  first  power 
of  the  frequency.  This  assumption  is  shown  to  be  valid  for 
at  least  a decade  of  frequency  by  the  attenuation  measure- 
ments made  in  the  large  tank.  The  dashed  line  in  Fig- 


ure VII-5  shows  how  well  such  a determination  of  the  at- 
tenuation fits  the  amplitude  data  at  a single  frequency. 

V 

The  log  decrement  is  obtained  from  the  measured  change  in 
attenuation,  Aa,  by 


Aa  V 
Af  8.686 


VI 1-3 


where  V is  the  shear  wave  velocity  and  8.686  is  the  conver- 
sion factor  between  logg  and  logig,  that  is,  a(dB/m)  = 

8.686  a(nepers/m). 

Shear  wave  speed  and  attenuation  were  obtained  as 
functions  of  frequency  for  two  of  the  sediments,  both  when 
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dry  and  completely  saturated.  Neither  the  PC  sand  nor  the 
MH  "beads  exhibited  dispersion  over  the  total  range  of  ap- 
proximately .5  to  20  kHz.  Values  at  different  frequencies 
differed  roughly  by  as  much  as  5 m/sec,  but  did  not  display 
a repeatable  pattern. 

The  two  receivers  used  simultaneously  in  the 
large  tank  would  allow  the  attenuation  to  be  calculated 
without  further  disturbance  of  the  sediment,  provided  the 
responses  of  the  transducers  were  known.  Since  no  standard 
unconsolidated  medium  of  known  shear  wave  attenuation  was 
available,  and  since  the  transducer  sensitivity  is  depen- 
dent upon  the  wavelength  in  the  medium,  calibration  in  the 
sediment  under  investigation  was  required. 

A successful  calibration  was  obtained  for  the  dry 

V 

PC  sand  by  placement  of  the  two  receivers  at  equal  dis- 
tances from  the  projector.  At  equal  separations,  the  ener- 
gy loss  due  to  geometrical  spreading  and  sediment  attenua- 
tion should  be  the  same  for  each  transducer.  The  measured 
amplitude  difference  therefore  represents  the  difference  in 
transducer  response.  An  average  correction  factor  was  ob- 
tained from  readings  at  several  equally  spaced  settings  and 
at  several  frequencies.  The  transducers  were  then  placed 
at  different  separations.  The  measured  attenuations. 
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corrected  for  the  transducer  response  and  assumed  spherical 
spreading,  are  shown  in  Figure  VII-6.  The  bars  on  the 
graph  are  one  standard  deviation  of  the  calibration  mea- 
surements divided  by  the  distance  between  the  transducers. 
The  open  points  were  obtained  at  different  relative  separa- 
tions than  the  closed  points.  Both  the  large  and  small 
transducer  elements  yield  similar  results. 

A consistent  calibration  was  not  obtained  in  the 
saturated  PC  sand  nor  the  MH  beads.  The  relative  amplitude 
difference  between  the  two  transducers  depended  on  both  the 
placement  and  the  frequency  and  was  not  closely  repeatable. 
Data  for  a single  transducer  showed  the  same  scatter  as  in 
Figure  VII-5.  It  was  again  necessary  to  introduce  the  con- 
cept of  a variable  transducer  coupling  in  order  to  extract 
meaningful  estimates  of  the  attenuation. 

With  the  amplitude  of  a single  transducer  ex- 
pressed by  Eq.  VII-2,  the  combination  of  two  receivers  at 
distances  xi  and  xj  gives 

AAj(xi)  - AA2(x2)  = ATi  - AT2  + Aa(xj-X2 ) , VII-U 

where  A represents  the  difference  between  the  two  frequen- 
cies. Reversal  of  the  separations  gives 


/ 


I 

I j 


13U 

AA 1 ( X2 ) - AA2 ( x 1 ) = ATj  - AT2  - Aa(xj-X2)  • VII-5 

Subtracting  these  two  expressions  yields  the  following 
equation  for  the  change  of  attenuation  between  the  two  fre- 
quencies considered. 

Act  = [ AAi  ( xj  ) -AA2  ( X2  ) -AAj  ( X2  )+AA2  ( xi  ) ^Ax  . VII-6 

The  variation  of  k * Ao/Af  indicates  the  frequen- 
cy dependence  of  the  attenuation.  A constant  k would  imply 
that  the  attenuation  was  proportional  to  the  first  power  of 
the  frequency.  For  the  dry  MH  beads,  a plot  of  k versus 
frequency  is  fairly  constant.  The  values  for  the  saturated 
sand  and  beads  are  scattered. 

If  a single  value  of  the  absolute  attenuation  was 
available,  then  the  changes  of  attenuation  calculated  from 
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Eq.  VII-6  would  yield  a plot  of  attenuation  versus  fre- 
quency similar  to  Figure  VII-6.  A reference  value  of  at- 
tenuation was  chosen  from  the  equation  a = kf  for  a fre- 
quency pair  whose  calculated  value  of  k matched  the  average 
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k for  all  such  pairs.  The  results  are  shown  in  Fig- 
ures VII-7*  VII-8,  and  VII-9.  Each  point  on  these  graphs 
was  determined  from  eight  separate  amplitude  measurements 
in  accordance  with  Eq.  VII-6.  Confidence  limits  for  the 

j 

curves  are  hard  to  assess,  particularly  in  terms  of 
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FIGURE  VI 1-9 

SHEAR  WAVE  ATTENUATION  versus  FREQUENCY 
SATURATED  MH  BEADS 
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absolute  values.  However,  the  results  were  repeatable  for 
both  sizes  of  transducer  elements  and  should  therefore  re- 
present a fair  picture  of  the  true  variation  of  attenuation 
with  frequency.  The  lines  with  unity  slope  depict  a linear 
dependence  of  attenuation  on  frequency.  They  are  included 


as  a reference  and  do  not  necessarily  represent  the  best 
possible  fit  to  the  data. 

A summary  of  the  acoustic  properties  of  the  sedi- 
ments as  determined  in  the  present  report  are  presented  in 
Table  II. 
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TABLE  II 


SEDIMENT  ACOUSTIC  PROPERTIES 
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CHAPTER  VIII 

DISCUSSION  OF  RESULTS 

The  compress' onal  wave  speeds  and  attenuations 
presented  in  Table  II  for  the  saturated  sediments  are  in 
accord  with  the  findings  of  other  investigators.  The  speed 
is  essentially  independent  of  grain  size  for  the  graded 
glass  heads.  Part  of  the  variation  that  does  exist  can  be 
attributed  to  changes  in  temperature  for  which  no  correc- 
tion has  been  applied.  The  average  value  agrees  with  a 
speed  of  1,850  m/sec  reported  by  Wylie  et  al . (1956)  for 
seven  sizes  of  glass  beads.  Speeds  observed  in  the  natural 
sands  are  slightly  lower  than  those  in  the  beads.  This  is 
expected  since  glass  has  a higher  bulk  modulus  than  quartz. 

The  damping  shows  a smooth  decrease  with  increasing  grain 
size  as  expected  for  sand  size  particles  (Hamilton,  1972). 

The  damping  is  smaller  in  the  angular  PC  sand  than  in  the 
other  samples. 

The  three  values  of  compressional  wave  speed  ob- 
tained in  the  dry  glass  beads  are  also  independent  of  grain 
size.  The  values  agree  with  longitudinal  speeds  of 
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98  to  113  m/sec  found  in  beads  by  Schmidt  (1951*).  Wylie 
et  al.  (1956)  estimated  a speed  of  UoOm/sec  in  dry  beads 
but  disclaim  reliance  on  this  figure  due  to  poor  arrivals 
and  scatter  in  repeated  measurements.  Such  a speed  corres- 


ponds to  values  found  in  the  present  investigation  whdn 
very  slight  frame  pressures  were  applied. 

Little  confidence  should  be  placed  in  the  quoted 
values  of  compressional  wave  damping  in  the  dry  materials 
due  to  the  difficulties  with  transducer  coupling  and  the 
resulting  poor  signal  quality.  The  values  do  agree  reason- 
ably well  with  decrements  obtained  by  Nyborg  et  al . (1950) 
of  roughly  1.2  to  1.6  and  by  Koltonski  and  Malecki  (1958) 
of  .23  to  .42  for  dry  sands  using  pulse  transmission  tech- 
niques. Recall  that  the  attenuations  observed  by  these  in- 
vestigators  did  not  vary  linearly  with  frequency.  Longi- 
tudinal decrements  obtained  in  dry  materials  at  small  con- 
fining pressures  are  considerably  lower.  For  example,  Pil- 
beam  and  Vaisnys  (1973)  quote  a Q value  that  corresponds  to 
a decrement  of  .07  for  longitudinal  vibrations  in  a sample 
of  MH  beads  at  .095  MPa  with  a speed  of  380  m/sec.  This 
value  did  not  vary  with  frequency. 


The  shear  wave  speeds  in 


Table  II  show  a greater 


percentage  variation  with  grain  size  than  do  the 
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compressional  wave  results.  However,  the  variation  fol- 
lows no  discernable  trend  with  either  the  grain  size  or  the 
porosity.  The  dry  sediment  speeds  follow  the  same  pattern 
hut  at  higher  values. 

The  shear  wave  decrements  for  the  saturated  heads 
decrease  slightly  with  increasing  grain  size.  Of  the  nine 
shear  wave  decrements  given  in  Table  II  for  the  saturated 
samples,  four  fall  within  the  range  of  .3  ± .15  suggested 
by  Hamilton  (1976b)  as  typical  for  surficial  sands.  The 
others  are  higher. 

A decrease  in  shear  wave  speed  with  saturation 
was  anticipated  from  the  theoretical  considerations  out- 
lined in  Chapter  III  and  from  the  experimental  results  for 
rocks  and  pressurized  sands  reviewed  in  Chapter  II.  All  of 
the  samples  studied  in  the  present  report  exhibit  this  be- 
havior. An  increase  in  shear  wave  attenuation  with  satura- 
tion was  also  expected.  This  was  observed  in  all  of  the 
natural  sands  and  in  the  MH  beads  in  the  large  tank.1 

1The  log  decrement  of  the  dry  PC  sand  in  th** 
small  tank  is  slightly  higher  than  that  in  the  saturated 
sample.  However,  the  attenuation,  that  is,  the  loss  in 
amplitude  with  distance  rather  than  wavelength,  is  less. 
Also,  difficulties  with  the  pulse  shape  were  encountered  in 
the  small  tank  at  high  velocities. 
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On  the  other  hand,  the  glass  beads  all  display  larger  at- 
tenuations when  dry  than  when  saturated.  This  difference 
is  apparently  due  to  the  instability  of  the  dry  glass  beads 
in  the  small  tank. 

The  extreme  sensitivity  in  the  small  tank  of  the 
dry  glass  bead  measurements  to  very  minor  physical  distur- 
bances makes  it  unlikely  that  the  measured  values  should  be 
taken  as  the  frame  parameters  when  the  same  beads  are  sat- 
urated. That  is,  the  contact  stresses  between  the  dry 
beads  in  the  wooden  box  probably  differ  substantially  from 
those  in  the  saturated  beads  in  the  small  aluminum  tank 
even  though  the  sediment  was  at  the  same  depth  in  each 
case.  Presumably,  the  increase  in  overburden  pressure  in 
the  large  tank  imparts  enough  rigidity  to  the  bead  packing 

V 

that  the  frame  parameters  in  that  instance  change  only 
slightly  with  saturation.  The  angular  sand  grains  are  less 
subject  to  changes  in  particle  arrangement  due  to  minor 
vibrations  and  are  also  less  susceptible  to  static  electric 
charges.  The  sand  frame  parameters  should  therefore  be 
less  influenced  by  saturation. 

While  the  exact  influence  on  the  frame  of  grain 


angularity  and  overburden  pressure  cannot  be  determined 
from  the  few  measurements  shown  in  Table  II,  important 
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Table  II  yield  k = .012  ± .005  for  the  beads  and  k = 

.05  ± .03  for  the  sands.  On  this  basis,  the  above  equation 
is  not  a particularly  accurate  means  of  estimating  shear 
wave  attenuation. 

A linear  increase  in  attenuation  with  freqeuncy 
has  been  found  in  the  large  tank  for  both  the  dry  MH  beads 
and  the  dry  PC  sand  over  approximately  a decade  of  frequen- 
cy (Figure  VII-6  and  Figure  VII-8).  A linear  increase  with 
frequency  is  also  compatible  with  the  results  obtained  in 
the  saturated  sediments,  although  the  scatter  in  the  mea- 
surements allows  other  interpretations.  In  particular,  it 
is  of  interest  to  compare  the  results  with  the  predictions 
of  the  viscous  loss  model  outlined  in  Chapter  III.  At 
first  glance,  the  observed  behavior  satisfies  the  major 

V 

theoretical  predictions.  That  is,  the  shear  wave  speeds 
decrease  and  the  attenuations  increase  with  water  satura- 
tion. Also,  the  attenuation  of  the  frame  can  be  character- 
ized by  a constant  log  decrement. 

Figures  VIII-1  and  VIII-2  compare  the  attenua- 
tions predicted  by  Eq.  III-19  with  the  experimental  data 
shown  in  earlier  figures.  The  frame  parameters,  y and  y', 
were  taken  from  the  dry  sediment  measurements.  The  pore 
size  parameter  necessary  for  the  determination  of  the 
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FIGURE  VI 1 1-1 

SHEAR  WAVE  ATTENUATION  COMPARED  WITH 
PREDICTED  VISCOUS  LOSS  - PC  SAND 
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FIGURE  VIII-2 

SHEAR  WAVE  ATTENUATION  COMPARED  WITH 
PREDICTED  VISCOUS  LOSS  - MH  BEADS 
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effective  density  was  obtained  from  Eq.  IV-4  and  the  data 
in  Table  I.  The  predicted  increase  in  attenuation  due  to 
viscous  loss  is  in  general  agreement  with  the  data  for  both 
the  MH  beads  and  the  PC  sand.  The  predicted  dispersion  in 
each  case  is  less  than  the  scatter  in  the  velocity  measure- 
ments. However,  the  predicted  speeds  are  slightly  higher 
than  those  actually  measured.  Minor  adjustments  in  the 
frame  parameters  can  result  in  better  agreement. 

The  general  agreement  of  the  measured  attenua- 
tions with  those  predicted  by  the  thee  ;tical  model  implies 
that  relative  movement  between  the  frame  and  pore  fluid  oc- 
curs during  an  acoustic  disturbance.  This  in  turn  suggests 
that  theories  which  assume  perfect  coupling  between  the 
frame  and  the  fluid,  such  as  that  of  Gassman  (1951),  may  be 
inadequate  for  a complete  description  of  sediment  acoustic 
properties.  This  contention  is  supported  by  the  shear  wave 
velocity  measurements  of  Hardin  and  Richart  (1963)  and 
Domenico  (1977)  who  found  that  the  saturated  speeds  in 
sands  were  greater  than  those  predicted  by  the  addition  of 
the  total  mass  of  the  water. 

If  the  viscous  loss  is  indeed  an  important  mech- 
anism, then  shear  wave  decrements  measured  at  a single 
frequency  cannot  be  confidently  extrapolated  to  frequencies 
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decades  removed  on  the  basis  of  a simple  linear  dependence 
of  attenuation  on  frequency.  Furthermore,  this  statement 
should  apply  equally  to  compressional  vave  measurements. 
This  in  turn  may  explain  the  variation  in  the  frequency  de- 
pendence of  compressional  vave  attenuation  reported  in  the 
literature  for  sands. 

The  relative  importance  of  the  viscous  loss  term 
depends  on  the  magnitude  of  the  frame  loss  and  the  ease 
with  which  the  fluid  moves  with  respect  to  the  frame.  The 
highly  permeable  sediments  studied  in  the  present  report 
were  specially  chosen  so  that  the  expected  viscous  loss 
would  be  comparable  to  dissipation  in  the  frame.  This  may 
seldom  be  the  case  in  natural  marine  sediments.  The  frame 
loss  is  expected  to  dominate  in  low  permeability  sediments 
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such  as  silts  and  clays.  Also,  increased  rigidity  due  to 
overburden  pressure  will  cause  both  the  frame  and  the  vis- 
cous losses  to  decrease,  although  not  necessarily  in  pro- 
portion. At  higher  pressures,  the  relative  displacement  of 
the  fluid  and  the  frame  may  be  negligible  even  in  sands.  A 
definite  statement  on  the  practical  importance  of  viscous 
loss  awaits  extended  measurements,  preferably  in  situ. 


CHAPTER  IX 

CONCLUSIONS 

The  data  contained  in  the  present  report  appar- 
ently represent  one  of  the  first  major  attempts  to  measure 
the  viscoelastic  shear  wave  properties  of  saturated,  uncon- 
solidated media  without  application  of  external  confining 
pressure.  Relevant  conclusions  can  therefore  he  drawn  not 
only  from  the  measurements  themselves  but  also  from  the 
measurement  techniques.  These  include  the  following. 

1.  High  quality  shear  wave  signals  can  be  gen- 
erated and  detected  by  bender  elements  in  low  rigidity  un- 
consolidated sediments.  The  broadband  characteristics  of 
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these  transducers  allow  speeds  and  attenuations  to  be  mea- 
sured over  more  than  a decade  of  frequency  with  the  same 
device.  However,  further  development  can  lead  to  addition- 
al refinement.  In  particular,  a better  understanding  is 
needed  of  the  transducer  response  to  a change  in  sediment 

1 

load  and  of  the  coupling  of  energy  into  the  sediment. 

2.  Control  of  minor  stresses  that  influence  the 
rigidity  of  the  sediment  is  necessary  for  a successful  lab- 
oratory program  of  shear  wave  measurement.  Both  the  finite 
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size  of  the  test  chamber  and  the  method  of  compaction  can 
influence  the  packing  of  sand  sized  particles.  Techniques 
which  exert  a controlled  external  frame  pressure  can  avoid 
some  of  the  problems  associated  with  minor  disturbances  of 
the  sediment.  However,  the  resulting  measurements  should 
not  be  considered  as  appropriate  to  surficial  sediments 
without  additional  confirmation. 

The  150  liter  sediment  tank  used  in  the  present 
study  represents  a convenient  size  for  laboratory  measure- 
ments in  sand.  There  is  sufficient  sediment  mass  to  reduce 
the  influence  of  minor  vibrations  on  the  packing,  and  yet 
the  mass  is  not  too  large  to  vibrate  at  amplitudes  nec- 
essary to  reach  minimum  porosity.  The  tank  is  also  small 
enough  to  allow  deaeration  of  the  entire  sediment  at  once. 

3.  An  irregular  porous  foam  liner  has  been  found 
to  be  an  efficient  absorber  of  shear  wave  energy.  This  al- 
lows interference  from  reflections  to  be  minimized. 

4.  Minor  changes  in  the  physical  properties  of  a 
sediment  exert  a larger  influence  on  shear  wave  propagation 
than  on  compressional  wave  parameters.  Shear  waves  there- 
fore offer  a greater  potential  for  acoustic  characteriza- 


tion of  sediments. 


5.  Increases  in  grain  angularity  and  overburden 
pressure  both  cause  significant  increases  in  sediment  ri- 
gidity and  decreases  in  shear  wave  damping. 

6.  Shear  wave  dispersion,  if  present,  is  small. 

7*  Shear  wave  attenuation  over  a limited  fre- 
quency range  is  adequately  described  by  a linear  variation 
with  frequency.  There  is  evidence,  however,  that  viscous 
interaction  can  be  a significant  loss  mechanism  in  highly 
permeable  sediments.  Theoretical  models  based  on  viscous 
interaction  predict  a change  in  log  decrement  with  frequen- 
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